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Abstract 

We give a review of the mathematical and physical properties of the cele- 
brated family of Calogero-like models and related spin chains. 
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1 Introduction 



The system of one dimensional particles with inverse-square pairwise interactions 
has assumed a "Jack-in-the-box" role in mathematical and theoretical physics over 
the last three decades: it pops up in various seemingly disparate situations, it is 
a recurring and unifying theme in high energy and condensed matter physics and 
constitutes the prime example of an integrable and solvable many-body system. Its 
range of applicability spans fluid mechanics, spin chains, gauge theory and string 
theory. This model has been labelled in the literature with various subsets and 
permutations of the last names of Francesco Calogero, Bill Sutherland and Jurg 
Moser. We shall simply call it the Calogero model, for brevity and in recognition of 
its original inventor. 

The Calogero model (and its many generalizations) has reached maturity, in the 
sense that its various aspects (classical, quantum, differential equation, operator, 
statistical, symmetry etc.) have been exhaustively analyzed and no new ground- 
breaking results seem to appear recently (I could be wrong!). It is unlikely, however, 
that it will be relegated to the shelves of mathematical physics for perpetuity. It has 
already enjoyed several revivals in its 35-odd year history and there are still open 
issues awaiting resolution. Chances are it has a few more surprises up its sleeve, 
their unveiling requiring, and offering back, new intuition. 

The purpose of this brief review is to expose the interesting physical and mathe- 
matical properties of the system and whet the appetite of readers for further study 
and research on the topic. It is not meant to be exhaustive, or rigorous, or all- 
encompassing: there are excellent extensive review articles that can serve this pur- 
pose. The hope is that the small size, physical slant and stress on newer develop- 
ments will make this narrative an accessible and motivational first piece. 

2 Basic features of the Calogero model 
2.1 Introducing the model 

The prototype of the model is the system of identical particles scattering on the line 
with inverse-square interaction potentials, as first introduced by Calogero |IJ. Its 
hamiltonian (in an obvious notation) is 



The particle masses m have been scaled to unity. 

The motivation for such a model comes from various perspectives. This is the 
only scale-free two-body potential that one can have: quantum mechanically, the 
potential scales like the kinetic term and therefore the coupling g becomes dimen- 
sionless (in T% — 1 units). The inverse-square potential would arise as a centrifugal 
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term in higher dimensions and is 'borderline' from the sense of stability: anything 
stronger than that at short distances and the particles would 'fall' into each other. 
Finally, from the long-distance point of view, the inverse square potential is border- 
line for statistical mechanics phase transitions: a stronger potential leads to phase 
transitions, while a weaker one does not. The Calogero model 'straddles the line' in 
many respects. 

To make the system bound, one has to introduce a kind of external 'box'. One 
way of doing this is to include an external harmonic oscillator potential that confines 
the system. This is nice, since it does not spoil the basic features of the system and 
leads to an integrable model, the harmonic Calogero model. Extracting thermody- 
namics from this system, however, is a bit tricky since the box is not homogeneous 
and the process requires careful scaling. Alternatively, one could put the system in 
a finite periodic box (whose length can be scaled to 2tt). The particles now inter- 
act through all the infinitely many periodic images of themselves and the two-body 
potential becomes 

V(X) = n^¥^^ = jd^f ^ 

This is the Sutherland model |2j. There are other versions of this class of models 
that we will not analyze. The original papers on the subject still consist essential 
reading [T]- [TTJ . The classic report [3] analyzes these systems in detail; the paper [TJ] 
describes relativistic analogs of these systems (known as the Ruijsenaars-Schneider 
model) that are not touched in the present review; the lectures cover many 
issues related to fractional statistics; and ^3] provides further mathematical results. 



2.2 Stability, hermiticity 

Classically the coupling constant g should be positive to ensure particles are not 
'sucked' into each other. Quantum mechanically the uncertainty principle works 
in our favor and the minimum allowed value for g is g = —\ (we put henceforth 
h — 1). This result can be derived by regularizing the singularity of the potential 
at Xi — Xj = and taking the limit, as done, e.g., in Landau and Lifshitz's classic 
quantum mechanics book. 

The hermiticity properties of the system for negative coupling constants has 
received a lot of attention and there are ways to make sense of values even less 
than g = —\- It is not our purpose to analyze these here, since it seems that the 
applications that arise naturally are the ones respecting the above condition. 

For later convenience, it is useful to parametrize g in the fashion 

g = £{£-!) (3) 

in which case the minimum value is naturally obtained for £ = |. Note, however, 
that there are generically two values of £ that give the same g, namely £ and £' = l—£. 
For — \ < g < 0, in particular, they are both positive. It is possible to argue that 
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only the value t > \ is relevant, although there are good reasons to retain both 
of them for the case of negative g. The possible use of that will become clear in 
the sequel. Note, further, that the above expression for g, when we reintroduce h, 
becomes g = £(£ — h). In the classical limit we have g = £ 2 , which can also be 
thought as the limit i >> 1. 



2.3 Properties of the classical motion 

The above system is classically integrable, which means that there are N integrals 
of motion in convolution, that is, N functions on phase space with vanishing Poisson 
brackets: 

{/ n ,/ m } = 0, n,m = l...N (4) 

For the scattering system (no external potential) I\ is the total momentum, J 2 is the 
total energy amd the higher I n are higher polynomials in the momenta also involving 
the two-body potentials. 

The integrals of motion derive from the so-called Lax matrix of the model. Con- 
sider the N x N hermitian matrices 

it N 1 1 

L jk =PjSjk + (1 - Sjk) — , A jk = £5 jk ^2 — + K$jk ~ 1)— (5) 

Xjk s=1 Xj s Xj k 

with £ 2 = g and Xjk = Xj — x^. Then, upon use of the equations of motion of the 
Calogero model, the evolution of the elements of the matrix L is 

L = i[L, A] (6) 

A pair of matrices satisfying this equation is called a Lax pair, with L the Lax 
matrix. What the above equation means is that the evolution of L is simply a 
unitary conjugation L — > ULU~ X generated by the (time-dependent) matrix A. So 
all the eigenvalues of L are conserved and the traces 

J„ = trL n , n = l,... TV (7) 

are constants of motion. Obviously l\ = J2iPi, while explicit calculation and the 
use of the identity 

i,j,k distinct l 3 %K 

(easily proved by multiplying numerator and denominator by Xjk and cyclically 
redefining the dummy indices i,j,k) shows that J 2 = 2H. I 3 and beyond generate 
new, nontrivial conserved quantities. 

It can further be shown that the above interals are in involution, that is, they 
have vanishing Poisson brackets. For the scattering system it is quite easy to give a 
physical proof: at t — > oo the particles fly far away from each other and L becomes 
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diag(pi), the off-diagonal elements asymptotically vanishing. The integrals of motion 
become simply 

In = £ K (9) 

i 

with ki the asymptotic momenta. Obviously these are in involution. Since the 
Poisson brackets of conserved quantities are also conserved, they must vanish at all 
times. 

The key interesting property of the above model, which sets it apart from other 
merely integrable models, is that, both classically and quantum mechanically, it 
mimics as closely as possible a system of free particles. We shall give here an 
overview of these properties, without insisting on proofs, and will come back to 
their derivation in subsequent sections. 

Let us first look at its classical behavior. The motion is a scattering event. 
Asymptotically, at times t = ±00, the particles are far away, the potentials drop 
off to zero and motion is free. When they come together, of course, they interact 
and steer away from their straight paths. Interestingly, however, when they are 
done interacting, they resume their previous free paths as if nothing happened. 
Not only are their asymptotic momenta the same as before scattering, but also the 
asymptotic positions (scattering parameters) are the same. There is no time delay 
of the particles at the scattering region. The only effect is an overall reshuffling of 
the particles. Thus, if one cannot tell particles from each other, and if one only 
looks at scattering properties, the system looks free! 

Similar behavior is exhibited by the harmonic Calogero and the Sutherland 
model. The motion of the harmonic model is periodic, with period determined by the 
harmonic oscillator potential. The particles seem to revolve around ghostly paths of 
would-be non-interacting particles in the external potential, their interaction creat- 
ing deviations but never completely upsetting the system. For the Sutherland model 
the particles seem bound to a set of free particle paths corresponding to a given set of 
momenta, deviating from them as they come near each other and resuming them as 
they separate. The momenta of the free paths (called 'pseudomomenta') are related 
to the integrals of motion of the system. 

Is there a concise manifestation of the coupling constant g in the classical motion 
of the particle, except its details? It turns out that there is a particularly neat result 
involving g: the classical action of two particles scattering off each other is the same 
as the action of free particles with the same asymptotic momenta, but diminished 
by 

AS = n^g = 7t£ (10) 

Semiclassically, the action corresponds to the phase of the wavefunction for the 
scattering process. This suggests that the scattering phase shift differs by a non- 
dynamical constant from the corresponding free particle result. As we shall see, this 
is fully born out in quantum mechanics. 
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2.4 Properties of the quantum system: fractional statistics 



The above behavior carries over to quantum mechanics. The asymptotic scattering 
momenta are the same before and after scattering. The classical fact that there 
is no time delay translates into the quantum fact that the scattering phase shift 
is independent of the momenta (remember that the time delay is the momentum 
derivative of the phase shift). Thus, it can only be a function of the coupling constant 
and the total number of particles. It is, actually, a very suggestive function: 

So the phase is simply £n times the total number of two-body exchanges that would 
occur in the scattering of free particles. This squares with the previous classical 
result: the classical limit corresponds to £ 1, so ^fg ~ £ and we get AS = 9 SC . 

We can interpret the above property as the fact that Calogero particles are 
essentially free but obey generalized statistics, as defined by their scattering phase 
shift [7 r Clearly the case £ = would correspond to free bosons and £ = 1 to free 
fermions (for these two values the potential vanishes and the system is, indeed, free). 
For any other value we can say that the system has fractional statistics of order £. 
Note that, in contrast to anyons, £ is not a phase: it can be fractional but also bigger 
than one ("superfermions"?). 

A word on the permutation properties of this system is in order. The inverse- 
square potential is quantum mechanically impenetrable, and thus the 'ordinary' 
statistics of the particles (symmetry of the wavefunction) is irrelevant: if the particle 
coordinates are in one of the N\ ordering sectors they will stay there for ever. The 
wavefunction could be extended to the other sectors in a symmetric, antisymmetric 
or any other way, but this is irrelevant for physics. No interference between the 
sectors will ever take place. All states have a trivial N\ degeneracy. (That is to 
say all states in all irreps of Sjv have the same physical properties.) Permutation 
statistics are therefore irrelevant and we can safely talk about the effective statistics 
as produced by their coupling constant £. 

Let us also review the properties of the confined systems. In the presence of an 
external harmonic potential of the form 

V = Y\^ (12) 

i 

the energy spectrum of a system of uncoupled particles would be 

E = -co + riiU (13) 
1 i 

The Hi are nonnegative integers satisfying 

Tit < ... < njv for bosons (14) 
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rii < . . . < tin for fermions (15) 
The actual spectrum of this model is 

E = yu; + £ N{N ~ 1] uj + riico (16) 

i 

with rii being 'excitation numbers' obeying bosonic selection rules: n« < Ui+i- Defin- 
ing the 'pseudo-excitation numbers' 

Ui = Hi + (i - l)t (17) 

we can check that the expression of the spectrum in terms of the fii is identical to 
the free one (|1H|) but with the quantum numbers now obeying the selection rule 

rk < - I (18) 

This is a sort of exclusion principle that requires the particle quantum numbers to 
be at least a distance i apart (as contrasted to 1 for fermions and for bosons). 
Again, a generalized statistics interpretation is manifest |17j . 

Let us clarify that the above numbers are no more integers. They do, however, 
increase in integer increments. The rule is that the ground state is determined by 
the minimal allowed nonnegative values for fii obeying (fT%J) while the excited states 
are obtained by all integer increments of these values that still obey (j!8j) . 

The periodic (Sutherland) model has similar properties. Its spectrum is 

E = J2 k + <E(*i - h) + e N{N l - l) (19) 

with the 'momenta' fcj being integers satisfying bosonic rules: ki < k i+ i. This 
looks rather different than the corresponding free expression (for i = 0). Defining, 
however, again 'pseudomomenta' 

Pi = ki + e(i- (20) 



2 

we can check that the expression for the spectrum becomes 

£ = E \v\ (21) 

i 

that is, the free expression. The pseudomomenta satisfy 

Pi < Pi+i - £ (22) 

that is, the same selection rule as the before. Again, we observe a generalization of 
the fermionic and bosonic selection rules corresponding to statistics £. The ground 
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state is the minimal (nearest to zero) numbers satisfying f!22l) while excitations 
correspond to integer increments (or decrements) of the pi still satisfying (|22|). 

Let us note that the above rule for £ = 1 reproduces the fermionic spectrum of 
particles with periodic boundary conditions for odd N and anta-periodic ones for 
even N: in the odd (even) N case the momenta are quantized to (half-) integers. 
This has a natural interpretation: when we take a particle around the circle it goes 
over N — 1 other particles. If we require the phase shift of the wavefunction in this 
process to agree with the minus signs picked up from the N — 1 fermion exchanges we 
recover the previous rule. We stress that, for free particles, this is not a consistency 
requirement but rather an aesthetic rule. At any rate, this is what the Sutherland 
model chooses to do! 

In conclusion we see that the Calogero model can be though of as a system of 
particles obeying generalized statistics. This manifests in terms of the scattering 
phases and, most significantly for statistical mechanics, through a peculiar 'level 
repulsion' of their quantum numbers generalizing the Fermi exclusion principle. 

2.5 Large-N properties of the Calogero model and duality 

Let us examine, now, the properties of the Calogero model as the number of particles 
grows large. At zero temperature, a non-interacting fermion system would form 
a Fermi sea. The corresponding 'Fermi surface' in one dimension degenerates to 
points. For the system in an external harmonic potential there is just one point 
corresponding to the highest excitation rip — N — 1. For the free periodic system 
we would have two Fermi momenta at ±pp = =t^r^ Excitations over this ground 
state are, then, conveniently classified in terms of particles (a filled Fermi sea with 
an isolated particle above or below) and holes (a filled sea with one unoccupied state 
inside it). 

Interestingly, the Calogero model presents a similar picture. The qualitative 
features of both the Calogero and the periodic Sutherland model are similar, so we 
pick the latter as most closely representing a gas of free particles in a box. From 
(J22|) above we see that the ground state also forms a 'pseudo-Fermi sea' (or should 
we call it a 'Luttinger sea'?) with Fermi levels rescaled by a factor £: p F = £^f^- 
Its minimal excitations are analogous to the ones of a Fermi sea, but not quite: 

• A particle would be an isolated occupied pseudomomentum above or below a 
completely filled sea of pseudomomenta. Particles 

-are excited in units of 1 (the increments of their pseudomomentum). 

-take up a space £ in pseudomomentum (since they cannot be 'packed' closer than 

£ units apart). 

• A hole would be an isolated empty space inside an otherwise occupied sea. 
Interestingly, the minimal such excitation is not obtained by removing one particle 
from the sea, but rather by incrementing all pseudomomenta of the sea above the 
place where we want to create the hole by one unit. Holes 
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-are excited in units of i. Indeed, since the distance of pseudomomenta in the sea 
is £, the possible positions of the hole are at distances I apart, 
-take up a unit space in pseudomomentum. Indeed, incrementing all pseudomo- 
menta above a given place in the sea by two units creates two holes in that place, 
and so on; by locally reshuffling pseudomomenta we can then separate these holes. 

Note that holes are not antiparticles. Removing a particle for the sea creates a 
gap of i spaces and, from above, I holes. The correspondence is 

1 particle ~ — t holes (23) 

We already observe a sort of duality between the two types of excitations. This can 
be summarized as 

particle <-> hole , I «-> — , p <-> £p (24) 

Under the above, the spectrum of excitations of the model remains invariant. This 
is the simplest manifestation of a coupling-constant duality that goes over to the 
correlation functions and Green's functions of the model [rH | ITH | l2T H l21j. Obviously, 
this duality is spoiled by nonperturbative effects, since holes are confined within the 
sea while there is no 'ceiling' for particle excitations. 



3 Particle symmetries and the Calogero model 

So far we have talked about the Calogero model and its properties considering it a 
'given' system, with only circumstantial motivation. It is now time to 'derive' the 
model, in the sense of obtaining it by starting from a set of principles. 

The angle we will take is the one of the description of a set of indistinguishable 
particles and their symmetries. Field theorists are used to the idea of viewing 
particles as field quanta, or as representations of the Poincare group. The many- 
body (first quantized) point of view, however, affords a perspective revealing the 
Calogero model as the natural generalization of free indistinguishable particles. 

The main idea is that identical particles admit the permutation group as a dy- 
namical symmetry commuting with the hamiltonian, while indistinguishable parti- 
cles elevate the permutation group to a (discrete) gauge symmetry. In plain words, 
configurations in which particles are permuted correspong to the same physical state 
and constitute 'gauge' copies of the system. 

There are two distinct ways to deal with a gauge system: 

• Reduce the system to a set of gauge invariant observables 

• Realize the gauge symmetry as a symmetry of the Hilbert space and impose gauge 
constraints on the states 

We shall explain below how each of these will lead to (versions of) the Calogero 
model. 



10 



3.1 Reducing to gauge invariant observables 

In the sense exposed above, phase space particle coordinates not physical 

observables since they insist on assigning a label to each particle; they are not 
permutation invariant. A set of invariants can be constructed in terms of symmetric 
functions of the above coordinates. Such a set is 

N 

In,m = E : *W : (25) 

8=1 

with n, m > and : ■ : denoting a specific ordering. E.g., the symmetric (Weyl) 
ordering between Xi and pi can be adopted, which also ensures the hermiticity of 

The above invariants are, in general, overcomplete. Even if we were to truncate 
the range of n, m to N, there would still be of the order of N 2 observables, while the 
number of independent basic operators is 2N. Classically, this reduncancy manifests 
in the existence of algebraic identities between various I m ,n- 

Quantum mechanically, overcompleteness translates into the presence of Casimirs 
in the algebra of I m ,n- Indeed, the l m ,n satisfy (a particular parametrization of) the 
so-called Wjy algebra. In the N — > oo limit their commutation relations can be 
conveniently repackaged into the 'sine algebra' ^3], by defining 

oo Lm„n 

E -rjlm,n (26) 

— 777 77 
m,n=0 " L - n - 

with k, q continuous 'Fourier' variables. Then, assuming Weyl ordering for the I m ,ni 
the I(k, q) satisfy 

[I(k, q), I(k', q')} = 2i sin W ~ k ' q i{k + k',q + q') (27) 

from which the commutators of I m n can be obtained by Taylor expanding I(k, q) 
and matching coefficients of k m q n . To lowest order in H we obtain 

[I m>n , I m , n ] = i(mn' - nm')/ m+m /_ lin+n /_i + 0(h 2 ) (28) 

which is the 'classical' algebra. (We have put h — 1, but lowest order in h 
corresponds to lowest order in m,n.) 

For finite N the algebra becomes nonlinear due to the presence of identities. 
Alternatively, we can keep the full tower of l m ,n arid effectively impose the identities 
as relations for the Casimirs of the algebra. Indeed, it is known that the above 
algebra admits a host of representations, corresponding to the underlying particles 
being bosons, fermions, parabosons or parafermions and various other possibilities. 

To see how the Calogero system emerges as one of these possibilities, concentrate, 
for the moment, to the special case of two particles ^E|- Their center of mass 
coordinate and momentum 

A = ^±^, P = Vi+V2 (29) 
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are certainly gauge invariant observables; they correspond to Ix,o and io,i- The 
relative coordinate and momentum, however, 



P1-P2 , qn v 
x = x x - x 2 , p = — - — (30) 



are not, since they are odd under permutation. We can form quadratic invariants 

as 

A 2 D Xp ~\- pX 2 / i\ 

A = x , B = , C = p (31) 

They correspond to 



A = 2/2,0 — 7 B = Ix,i — -(Ix,oIo,i + h,ih,o) , C = -J ,2 — (32) 

The relative variables A, B, C commute with the center of mass variables and close 
to the SL(2, R) algebra 

[A, B] = ilA , [B, C] = ilB , [A, C] = iAB (33) 

Classically, the above variables satisfy the constraint AC = B 2 . Quantum mechan- 
ically, this translates to the Casimir 

G = AC + CA _ b 2 (34) 

Physical Hilbert spaces correspond to irreducible representations of the algebra of 
A, B, C, along with X, P. 

The representation with vanishing Casimir corresponds to the original system of 
two bosons or two fermions, in which A, B, C can be realized as in ()31j1 . This is not, 
however, the only one. A nonzero value for G, corresponding to a quantum correction 
to the classical value, would also be a legitimate realization of the indistinguishable 
particle dynamics. Interestingly, unitarity mandates that G > — |, so we may 
parametrize 

G = £(£ - 1) (35) 

in direct analogy to the Calogero case. We observe that the representation with the 
above value of G can be realized as 

W, S = 5E±J^, c = P 2 + ^ (36) 

2 X z 

Effectively, the relative kinetic energy of the particles has acquired an inverse-square 
potential part. The free particle hamiltonian for this system would become 

rr P2 1 2 1 2 W ~ 1) 

This is the Calogero model! The Calogero coupling plays the role of the Casimir. 
The point is that the hamiltonian and other observables of the Calogero model 
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and the corresponding 'free' system (free of inverse-square interactions) in terms of 
physical variables X, P; A, B, C are identical. 

For more than one particles the approach can be extended and we again recover 
the Calogero model as one particular realization of the indistinguishabe particle 
system. 

The Hilbert space of the above realization, in terms of X{,pi, still provides a 
defining iV!-dimensional realization of the permutation algebra. For generic values 
of the Calogero coupling, however, this is inconsequential: the inverse-square poten- 
tial is impenetrable (both classically and quantum mechanically) and thus the N\ 
sectors do not mix and are physically equivalent. There is no need to decompose 
the defining representation into irreducible components. The 'normal' statistics of 
the particles (wavefunction symmetry properties) have become immaterial, being 
totally supplanted by the Calogero dynamics. Only for G — 0, that is, t = or 
1, the Calogero interaction vanishes and ordinary statistics come back into play. 
By analytic continuation of the wavefunction properties, we can map bosons to the 
I = case and fermions to the i = 1 case. 



3.2 Augmenting the symmetry: Matrix Model 

The opposite way to realizing a gauge system is to keep the original redundant for- 
mulation and impose gauge invariance as an operator relation on the states (Gauss' 
law). In terms of the single-particle phase space variables x = i,p = i the gauge sym- 
metry is the permutation (symmetric) group S = N. This leads to the well-known 
and analyzed cases of fermions, bosons and their parastatistics generalizations. 

A different approach is to start with an augmented system, in which both the 
dynamical variables and the gauge symmetry have been expanded, giving at the end 
the same gauge invariant degrees of freedom [22] • Specifically, we could formulate 
the particles in terms of the eigenvalues of an N x N matrix. There is no a priori 
ordering of these eigenvalues, so this certainly encodes identical particles. It is clear 
that the permutation symmetry of the problem has been promoted to the continuous 
symmetry of unitary conjugations of this matrix, which leaves the eigenvalues intact. 

This unitary conjugation symmetry is now the gauge group. States need not 
necessarily be singlets under this symmetry, however, just as in the previous sec- 
tion the Casimir G needed not vanish. We can simply select the Hilbert space to 
transform under an irreducible representation of the gauge group and identify states 
transforming to each other under the action of the group as a unique physical state, 
in direct analogy to the considerations that lead to parastatistics. 

This is the celebrated matrix model formulation, to be fully analyzed in the 
subsequent sections. It has found various applications in physics, the most directly 
related to identical particles, perhaps, being the noncommutative Chern-Simons 
description of the finite quantum Hall droplet [28]-|27j. 
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4 The hermitian matrix model 



4.1 Classical analysis 

Let us first examine a matrix model that parallels as closely as possible particle 
mechanics on the infinite line. The kinematical variable is a hermitian N x N 
matrix M and the lagrangian reads 

£ = tr|^M 2 -\/(M)| (38) 

V(x) is a scalar potential evaluated for the matrix variable M. 

Clearly the above has a time-translation invariance which leads to the conserved 
energy 

H = tr j^M 2 + y(M)J (39) 

Moreover, the action is invariant under time-independent unitary conjugations of 
the matrix M: 

M -> UMU' 1 (40) 
This nonabelian SU (N) symmetry leads to the conserved hermitian traceless matrix 

J = i[M,M] (41) 

where [ , ] denotes ordinary matrix commutator. These are the 'gauge charges' 
that, when fixed, will determine the particular realization ('statistics') of the indis- 
tinguishable particle system. But let us further analyze the implications of fixing 
these charges classically. 

We are interested in the dynamics of the eigenvalue of M, so we parametrize it 

as 

M = UAU- 1 (42) 

where U(t) is the unitary 'angular' part of the matrix and A(t) = diag{x\, . . .x^} 
are the eigenvalues. Clearly the conserved quantity J has to do with invariance 
under 'rotations' of the angular part of M and thus corresponds to the 'angular 
momentum' of U(t). We define the 'gauge potential' 

A = -U- l U (43) 

M, J and the lagrangian C become, in this parametrization, 

M = U (A + [A, A}) U- 1 (44) 
J = iU ([A, [A, A]]) U~ x = UKU -1 (45) 

C = tr{iA 2 + [A,A] 2 -\/(A)} (46) 

j N i N 

= -j^Z^i ~ 2^Zi x i — Xj) AijAji (47) 

i=l i,j=l 
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The matrix elements of A and K are related 



K jk = i [A, [A, A]] jk = ifa - x h yA jk (48) 
Finally, solving (|4*Hj) for Aj k and putting into (}4Tj) we obtain 

The first two terms are kinetic, coming from M 2 , while the last one is potential. 
Therefore, the hamiltonian H is 

Note that the eigenvalues are kinematically coupled by an inverse-square type poten- 
tial with the angular momentum degrees of freedom. The connection of the matrix 
model to the Calogero model along the lines presented here and below was first 
established in Also, the hamiltonian (J5U)l has been proposed independently of 
the matrix model as an S{7(A/")-generalization of the classical Calogero system j^J. 
We can now examine special cases: 

a) The most 'gauge invariant' sector is, of course, the one in which the angular 
momentum charges vanish, that is, J = 0. In that case, (|5T)|) for K = becomes the 
hamiltonian of non-interacting particles in an external potential V(x). This would 
be the case of 'standard' particles. 

b) For the next simplest case the angular momentum J should be as simple 
as possible without vanishing. Only its eigenvalues are really relevant, since we 
can always perform a time-independent unitary transformation V which would shift 
U — > VU and would rotate J — > VJV~ l . The simplest choice would be to take the 
eigenvalues of J to be equal. Unfortunately, this is not possible since the traceless 
condition would make them vanish. The simplest possible choice is to take all the 
eigenvalues equal to C except one, which would cancel the trace by being (1 — N)£. 
This can be written in terms of an arbitrary column A^-vector v as 

J = £(W - 1) , v ] v = N (51) 

in which case K becomes 

K = £(uu^ - 1) , u = XJ- x v (52) 
From we see that Ka = (no sum on i) and thus 

UiU* = 1 (no sum) (53) 
So the coefficient of the inverse-square potential in (joUj) becomes 

KijKji = lmu* iujU* = £ 2 (i ^ j) (54) 
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Finally, (JHOj) becomes 



1 



H = £ + E 7^— ^ + E ^) (55) 

This is the Calogero model! The potential strength g = i 2 is related to the conserved 
charge t. Quantum mechanically, picking this charge will amount to a choice of 
statistics. We also get, at this stage, an arbitrary external potential V(x). 

c) More general choices of J amount to more variety in its eigenvalues. KijKji 
now, in general, becomes time- dependent and the dynamics more complicated. We 
postpone the discussion for the quantum case where it will be shown that this 
corresponds to Calogero particles having also internal degrees of freedom. This 
will be a generalization of the discussion of the first section, with irreps of SU(N) 
substituting the irreps of Sn- 

Now that we have this new approach we can use matrix technology to demon- 
strate the integrability of the Calogero model [HE]. For V(x) = the matrix motion 
becomes free and M is conserved. The conjugation-invariant quantities 

I n = trM n (56) 

are also conserved and in involution (the matrix elements of M are momenta and 
have vanishing Poisson brackets). From (|44jl and ()48j) we have 

(UMU- 1 )^ = 5 jk x 3 -{l-5 3k )^^- (57) 

111 '?/ 

= 5 ]k x 3 -{l-5 3k )-^- (58) 

Xj x k 

(note the similarity of the above expression for UMU~ l with the Lax matrix defind 
in section 2.3). When the above expression is inserted in the trace I n = trM n clearly 
U drops and products of the form u(a* UjU* k . . . u* will appear which reduce to powers 
of £. Therefore, the I n reduce to expressions involving only Xj, Xj and the coupling 
constant £. These are the conserved integrals of the Calogero model. 

Starting from the matrix model the actual motion of the Calogero model can be 
obtained. The solution for M is 

M = B + Ct (59) 

for arbitrary matrices B, C. The conserved charge becomes 

J = i[M, M] = i[B, C] = il{uv) - 1) (60) 

By unitary transformations we can choose the phases of u such that Ui—1; choices 
for B, C, then, that satisfy (J60|) are 

U 

B jk = S jk q 3 - , C jk = 5 jk pj - (1 - 5 jk ) (61) 

Qj - Qk 
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qi and pi are the initial conditions for x% and X{ at time t = 0. Diagonalizing, then, 
M = B + Ct for the above B, C produces the motion of the system. Another choice, 
in which C is diagonal, is 

it 

Bjk = &jk a j + (1 _ fijkjT r~ j = Sjkkj (62) 

% "-ft 

fci and Oj are asymptotic momenta and impact parameters. For t — ► ±oo the off- 
diagonal elements of 5 produce a perturbation of order t -1 to the eigenvalues, so 
the motion is determined by the diagonal elements + kit alone. We recover the 
result that the motion at asymptotic regions is the same as if the particles were free. 

We conclude by proving that the matrix model is also integrable and solvable in 
the presence of a harmonic oscillator potential V{x) = hjj 2 x 2 . The non-hermitian 
matrix Q = M + iuoM evolves as 

Q(t) = e^Q(O) (63) 

and the matrix Q^Q is conserved. We leave it as an exercise to derive the conserved 
integrals and the motion of the corresponding Calogero problem. 

External potentials with up to quartic dependence on x also lead to integrable, 
although not so solvable, models [2H]- It is an open question to prove this is all there 
is, or to find yet more integrable potentials. 

Finally, we may wonder what restricts us to one dimension. We chose a model 
with one matrix, and its eigenvalues corresponded to coordinates of particles on 
the line. We could, indeed, start with an appropriate model with many matrices, 
which would reproduce particle motion in higher dimensions (2H1- The integrability 
and solvability properties of such extended models, however, are much less pleasant. 
The question of whether they represent a workable extension of identical particles 
remains open. 

4.2 The hermitian matrix model: quantum 

We will, now, perform the quantization of this system. We will do it first for the 
hermitian matrix model and subsequently for the unitary one. Each model has its 
own advantages and appears in different situations. 

Consider the hermitian model with a quadratic potential, 

C = tr (-M 2 - u 2 M 2 

with uj a scalar frequency. The above can also be written in terms of matrix elements 

£ = 1 E\\Mj*\ 2 -w 2 \M jk \ 2 (65) 
jk 

This is nothing but iV 2 harmonic oscillators. So the system is in principle trivial 
and solvable. All its nontrivial features emerge from the reduction to a subspace 
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corresponding to a fixed value for the 'angular momentum' J. Reducing M to its 
matrix elements is not beneficial for this purpose; we need to treat it as a matrix. 
We begin by defining a canonical momentum matrix conjugate to the 'coordinate' 

M 

P=^L = M (66) 
In terms of M and P the hamiltonian of the model becomes 

H = ti(-P 2 + -cu 2 M 2 ) (67) 



2 2 
The Poisson brackets are 

{M jk , P lm } = S jm 5 lk (68) 

Upon quantization, the matrix elements of M and P become operators and the 
above Poisson brackets become quantum mechanical commutators (not to be con- 
fused with matrix commutators). The commutator or Mj k and Pi m can be written 
conveniently by thinking of indices j, k as acting on a linear space 1 and indices I, m 
as acting on linear space 2. The 4-index symbol 

8jm&lk = (Ti 2 )jk;lm (69) 

acts on both spaces and is, in fact, the operator exchanging the two spaces. Denoting 
with X 1 and X 2 any matrix X acting on space 1 or space 2, respectively, T 12 satisfies 

X\T\2 = T12X2 , X 2 Ti 2 = T\2X\ (70) 

We also note the partial trace and unimodularity properties 

t ri T 12 = I 2 , tr 2 T 12 = h , T 2 2 = I 12 (71) 

with Ii, I2 and ii 2 denoting the identity matrix in space 1, space 2, or space 1 x 
2 respectively. With the above notation the quantum commutator of M and P 
becomes 

[M u P 2 ]=iT 12 (72) 

The complex matrix Q = M + iuM introduced before, and its conjugate Q\ are 
the matrix analogs of creation and annihilation operators. Defining matrix operators 
with the standard quantum normalization 

At = J=(P + iwM) , A=-±=(P-iuM) (73) 
\/2uj \/2uj 



they satisfy 

[A,4]=T 12 , [A u A 2 ] = [A\, A 2 ] = (74) 
By analogy with the standard harmonic oscillator, let us define the matrices 

L = A^A , R' = -AA^ (75) 
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The matrix operator L = A^A is normal ordered, since all (quantum) creation 
operators are to the left of annihilation operators. R', on the other hand, is not. 
We can define a new (quantum) normal ordered operator R, in which the matrix 
multiplication is first performed in the order AA^ and then the creation matrix 
elements are moved to the left. Specifically 

R jk = - : A ]S {A^) sk := {A r ) sk A JS (76) 

Using the commutation relations of A and A^ we see that R and R' are trivially 
related: 

R = R' + N (77) 
In terms of L and R, the hamiltonian can be written 

1 + + N 2 N 2 

H = -utr(AA f + A' A) = utrL + —to = -uotiR + — uj (78) 

From the basic commutators (|74|) we can infer the commutation relations: 

4] = T 12 A\ , [R u 4] = 4T 12 (79) 

as well as their hermitian conjugates. Tracing the first equation above with respect 
to space 1, we obtain 

[H, A j ] = uA ] (80) 

This means that any matrix element of A is a creation operator, creating one quan- 
tum of energy uj. This is hardly surprising, since the matrix model is, indeed, iV 2 
harmonic oscillators with identical frequency uj. 

The commutators of the matrix elements of L and R are calculated to be 

[LxM = {Lx-L 2 )T 12 
[Rx : R 2 \ = (Ri — R 2 )T\ 2 

[Li,R 2 ] = (81) 

The above is nothing but two commuting copies of the U(N) algebra in disguise. 
To see this, define the fundamental U(N) generators T a , a = 0, 1, . . . iV 2 — 1, with 
T° = I/y/N satisfying the normalization condition 

ti(T a T b ) = 5 ab (82) 

and the U(N) commutation relations 

[T a ,T 6 ] matrix = ^ abc T c (83) 

where [., .] matrix is a matrix commutator. Then expand L and R in the complete 
basis T A of hermintian N x N matrices: 

L a = tr{T a L) , L = T a L a (84) 

a 
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and similarly for R a . The (scalar) expansion coefficients L a , R a , upon use of (|81j) . 
satisfy 



[R a } R b ] 
[L a ,R b \ 



if abc L c 
if abc R c 




J5) 



Since L and R mutually commute, their powers and traces commute as well, so 

[tiL n ,trR m } = (86) 

On the other hand, using the fundamental commutation relation for A and A' ([740 
and the properties of T i2 f)70|71j) . we can relate the traces of L and R as 



triL? = tri^AO" = tr x [Ajtr 2 Ti 2 Ai(44 



n-l 



tratr 2 [4Ti 2 ^ a (^Ui 



>n— 1 



tr 



12 



,n-l- 



n-l 



$> 12 r 12 (A 1 4) fc T 12 (A}A 1 ) n - 1 - 



+ tr 



12 



Ti 2 (A 1 A\) n ~ 1 A 1 Al 



k=0 

n-l 



- X] tT(AA r ) k tr(A^A) n - 1 ~ k + tr(AA f ) n 



fc=0 



^(-l) fc tr(i? - N) k -HiL n - k + (-l) n tr(R - N) n 



57) 



fc=i 



Working recursively with the above relation, we can express traces of L entirely in 
terms of traces of R and vice versa; mutual commutativity of the two sets, therefore, 
translates into commutativity of the elements of each set, which shows that the 
quantities 

I n = tr(y4y4 f ) n (88) 

are conserved and in involution. 

The above proof, involving two commuting matrix operators L and M that are 
subsequently related, may seem a bit 'too slick'. Nevertheless, the involution of 
I n should be clear from ([81 jl : the I n are nothing but the Casimirs of the U(N) 
algebra generated by L a . I± = trL ~ H is simply the U(l) charge, while the rest 
are (related to) the higher SU(N) Casimirs. Since SU(N) is of rank N — 1 and 
has N — 1 independent Casimirs, we recover N commuting charges altogether. The 
higher charges I n , n > N, are related to the fundamental charges 7i,.../jv via 
nontrivial, nonlinear, iV-dependent relations. 



4.3 The hermitian matrix model: reduction and spectrum 

It remains to do the reduction to the 'gauge' sectors, that is, fix the charge J = 
i[M, M] ma trix and work out its implications for the quantum states of the system. 
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The 'angular momentum' J classically generates unitary conjugation of the ma- 
trix M; this implies that quantum mechanically it will become an SU (N) algebra. 
This can be seen explicitly: from its definition J can be expressed in terms of the 
operators defined in the previous section 

J = i:[M, P] matrix : = : [A\ A] matrix : = A* A - : AA* : = L + R (89) 

L and R are two commuting U(N) algebras, therefore their sum is another U(N) 
algebra. 

Classically J is traceless and therefore its U(l) part vanishes, making it an 
SU (N) matrix. Quantum mechanically this should still be true, since it generates 
the transformation 

M -> UMU- 1 (90) 

which has trivial U(l) part. This is ensured by the normal ordering of the above 
expreession for J. L and R satisfy trL = — tri? and so tr J = 0. 

Reduction of the system to particular values of J corresponds, quantum me- 
chanically, to fixing the representation of the SU (N) algebra J. So the system 
decomposes into sectors labelled by the allowed irreducible representations (irreps) 
of J. Further, states within each sector related via the action of J are identified as 
a unique physical state, since J is a 'gauge' symmetry. 

To identify these sectors we first need to identify the possible irreps for L and 
R. This can be done by examining their form. Take L, at first: it is nothing but 
N copies of the Jordan- Wigner bosonic oscillator construction of the U (N) algebra 
over the fundamental representation. This may require some explaining. 

Let R^p be any d-dimensional representation of a Lie group, where a, (5 — 1, . . . d 
label its matrix elements and a labels its generators. Define a set of d commuting 
creation and annihilation operators a a ,a) a 

[a a ,a P ] = 5 a/3 (91) 

The operator obtained by 'sandwiching' the matrix R A between the vectors a* and 
a 

G a = alR^ap (92) 

satisfies the commutation relations of the Lie algebra. Therefore, it provides repre- 
sentations of the algebra, imbedded in the Fock space of the oscillators. Specifically, 
it provides all representations generated by the fully symmetrized direct product of 
any number of representations R. This includes the singlet (the Fock ground state), 
R itself (the set of d states with excitation number one) etc. 

For the specific case of the U (N) Lie algebra with R the fundamental represen- 
tation F, we need iV oscillators and we get 

G a = a\T a jk a k (93) 
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Transforming G from the generator basis G a to the matrix basis G = J2 a G a T a , G 
simply becomes 

G jk = a)a k (94) 
In view of the above, the matrix elements of L can be written 

L jk = (A%A sk = A\jA sk (95) 

It is clear that j and k play the role of the fundamental indices of U (N), while s is a 
summation index that runs over N values; for each fixed s, A\j and A sk play the role 

of Oj and a k , respectively, while for different s all operators commute. So the above 
L is the direct sum of N independent (commuting) Jordan- Wigner realizations of 
U(N) over the fundamental. 

For each fixed s this realization includes, as explained earlier, all the totally 
symmetrix tensor products of the fundamental representation of U(N); that is, 
all irreps with a single row in their Young tableau. The direct sum of N such 
representations, however, includes irreps with up to N rows, which is the general 
case. We conclude that the spectrum of L spans the full set of irreps of U(N). The 
£7(1) charge (trL) is simply the total excitation number and is given by the number 
of boxes of the irrep. 

A Jordan- Wigner construction based on the antifundamental representation T a = 
— (T a )* = —(T a Y would lead to an expression in terms of N oscillators 

G jk = -a\ aj (96) 

Comparing with the expression for R 

R jk = -A\ s A js (97) 

we conclude that R is the direct sum of N independent Jordan- Wigner realizations 
of U (N) over the antifundamental, which again spans the full set of irreps of SU (TV). 
The U(l) charge, now, is the negative of the total excitation number. 

Representations of L and R are, in fact, constrained to be conjugate to each 
other. This arises due to their construction in terms of the same bosonic creation 
and annihilation operators, and is also manifest by the relation between their traces, 
which makes all their even (odd) Casimirs equal to (minus) each other. So the 
representations that J = L + R can carry are of the form r x f. Such representations 
include the singlet, the adjoint etc. and they have as common property that their 
Z N charge vanishes. We conclude that the irreps of J must have a number of boxes 
in their Young tableau that is an integer multiple of N. 

The above result for J could have been obtained by noticing that J itself is 
a bosonic Jordan- Wigner construction over the adjoint representation of SU(N). 
Clearly the adjoint and all its symmetric products have a number of boxes that is a 
multiple of N. 
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The quantum states of the model can be constructed by starting with the ground 
state |0 >, annihilated by all operators Aj k 

A jk \0 >= (98) 

and acting with any number of creation operators: 

h;3 2 , k 2 ; . . .j n k n >= A] iki A] 2k2 . . . At n JO > (99) 

For any c— number (classical) matrix $, the commutation relations 

[tr($L),A f ] = $A f 
[tr ($/?), A+j = -A f $ 

[tr($J),A f ] = $A^-A^ (100) 

imply that the first (left) index of the matrix (A')jk transforms in the fundamental 
(F) under L, the second (right) index transforms in the antifundamental (F) under 
R and that A* altogether transforms in the adjoint under J. The general state 
transforms in a representation of the form FxFx...xFxFx... with n Fs and 
n Fs appearing, corresponding the the number of free left and right indices in the 
state. 

Physical states are chosen by imposing the constraint that J is in a fixed represen- 
tation of SU (N), say r. To do this, we start from a generic state \j±, ki, j 2 , k 2 ; . . . j n k n > 
and contract the free indices with Clebsh-Gordan coefficients that project it to this 
representation. All states within the same representation r are gauge-equivalent and 
therefore represent a unique physical state. States corresponding to inequivalent 
copies of r, however, contained in the representation of the initial state correspod 
to distinct physical states. 

The above analysis can be done explicitly in the case of the representation r 
corresponding to the standard (spinless) Calogero model. Classically the matrix 
commutator i[M, M] = J is of the form £(vv^ — 1), where v is a vector of length 
squared equal to N. The representation rp corresponding to this is the fully sym- 
metric one with a number of Young tableau boxes equal to £N. 

The easiest way to realize this is by quantizing the classical expression J = 
£vv^ — i = ipip^ — t promoting the vector components ipjl^Vj to yet another set of 
harmonic oscillator creation operators 

[Ml] = s jh (ioi) 

In this way the operator 

Wik = -^j (102) 

realizes the U(N) algebra, again in a Jordan- Wigner construction over the anti- 
fundamental. Its representation content includes all fully symmetric products of 
the antifundamental, each of them represented in the subspace of fixed excitation 
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number n = ipjipj. The restriction of J to the symmetric representation re can be 
expressed by the condition 

J + + £ = o (103) 

This means that the representations carried by J and must be conjugate to 
each other so that their sum (composition) contain the identity, which fixes the 
representation of J to be a totally symmetric one. The c-number term i serves the 
purpose of subtracting the trace of J^, which is required by the tracelessness of J. 
It also fixes the U(l) part of J^, since tracing the above relation gives 

trJ^ + £N = - tptyj +£N = (104) 
j 

Wde recover the condition that n = £N is the total number of boxes in the Young 
tableau of J$ and thus also of J, fully fixing the desired representation. 

As we stated earlier, representations of J must have a number of Young tableau 
boxes that is an integer multiple of N. An important corollary of the above analysis, 
therefore, is that i must be quantized to an integer. This is a new feature of the 
Calogero system as deriving from the matrix model. It is not a general requirement 
for the quantum Calogero model, which is perfectly well-defined for fractional values 
of £. The source of this quantization is the enlargement of the symmetry group of 
the system from Sn to SU(N) for the matrix model: the enlarged symmetry has 
a global 'anomaly' which requires the quantization of the coupling constant £, in 
analogy with similar effects in gauge theory. 

Finally, states of the theory can now be constructed in terms of the vaccum 
state annihilated by A and if) upon the action of if)* and A*. The full state must 
be a singlet under J + J^, which means that all indices, including those of if)*, 
must be contracted. Gauge invariant contraction of indices can be done by matrix 
multiplication, but tracing or by multiplication of the vector ip* and the matrix A*. 

Each if>* in the state, therefore, will leave one index hanging, since there is no 
way to contract it with another if>* and multiplication with A^ still leaves one un- 
contracted index. The only way to contract these indices is using the only invariant 
tensor of SU(N) that includes all fundamental indices, namely the iV-fold antisym- 
metric tensor £j 1 ...j N - This means that -?/^s must come in multiples of N, recovering 
once more the condition that the total excitation number n = £N must be a mul- 
tiple of N. Further, operators contracted with e must all be distinct, otherwise the 
product would vanish due to antisymmetry. 

The generic form of the physical states is 



t 

> 



[trAtp [tr(At)f 12 . . . [tx{AYr h,., N ^W^)n ■ ■ ■ (^W"V 

(105) 

Higher traces of A* can be related to the first N ones, and other, more general ways 
of contracting operators with e tensors can be reduced to linear combinations of the 
above states. 
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This is an eigenstate of the hamiltonian H = uoti{A^ A) +ujN 2 /2 with its energy 
given by the total number of oscillators appearing in the state plus a constant. 
The -^-dependent prefactor contributes an energy equal to uo£N{N — l)/2. We 
conclude that the energy spectrum is 

E = u (X; m k + i^- 2 1 + — J (106) 

We recognize the excitation energies as those of N noninteraction bosons in a har- 
monic oscillator potential, expressed in terms of their collective excitations {my. 
represents the energy gap between the top k bosons and the next lower one on the 
single-particle oscillator spectrum). By standard bosonization arguments, the same 
excitation spectrum applies to a set of N noninteracting fermions in a harmonic 
oscillator (m& + 1) represents the energy gap between the top k fermions and the 
next lower one on the single-particle oscillator spectrum). So we may rewrite the 
above spectrum as 

E = a (jt* + l^ + %)=»(£*, + i) d07) 

where ri\ < HI2 < ...tin are single-particle bosonic excitation numbers while the 
'pseudoexcitation' numbers n« have been defined as 

fit = ru + (e + - 1) (108) 

We recover the spectrum of the Calogero model as exposed in a previous section, 
with the extra shift £ — > £ + 1. This is a quantum shift of the 'bare' parameter 
I appearing in the classical matrix model to the renormalized value I + 1 entering 
the Calogero model. In particular, the single sector I = of the matrix model 
corresponds to £ = 1 in the Calogero model, that is, fermions. The 'fermionization' 
of the eigenvalues of the matrix model due to the quantuym mechanical measure 
arising out of integratig out the angular variables of the matrix is a well-known 
effect. 

Other, more general representations of J can be dealt with in similar ways. As 
shall be explaind later, they correspond to Calogero particles with internal degrees of 
freedom ('spin'). We shall examine these cases in the context of the unitary matrix 
model. We have analyzed quite enough already the hermitian matrix model and we 
should leave something interesting for the unitary case! 

5 The unitary matrix model 
5.1 Classical analysis 

The hermitian matrix model works well for particles on the line but has trouble 
representing particles on periodic spaces. The most natural candidate for such 
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models would be a unitary N xN matrix U. Its eigenvalues are phases and naturally 
live on the circle. We start, therefore, with a lagrangian that represents the invariant 
kinetic energy on the space of such matrices: 

C = -UriU^Uf (109) 

A potential could in principle be included but we are interested in the translationally 
invariant case and will omit it. The treatment is similar as before, and we just 
summarize the relevant facts. 

The lagrangian is, in fact, invariant under separate left- and right-multiplications 
of U by time-independent unitary matrices and and there are two corresponding 
conserved matrix angular momenta L and R: 

U^VU: L = iUU- 1 (110) 
U^UW~ l : R = -iU^U (111) 

The unitary conjugation that preserves the eigenvalues corresponds to W = V and 
its generator is 

J = L + R = i[U,U~ 1 ] (112) 
The rest of the discussion is as previously. Parametrizing 

U = VAV- 1 with A = diag{e ix \ . . . e iXN } (113) 

the hamiltonian becomes, after a few steps, 

i z z i^j 4 sin 

where, as before, 

K = V- l JV (115) 

Choosing J = K = reproduces free particles on the circle, while choosing J = 
£(uu^ — 1) we obtain K^K^ = £ 2 and we recover the Sutherland inverse-sine-square 
model 

11 f 2 

H = y-x 2 t + -y j — —— (lie) 

^2 1 2^4sin 2 ^ 1 ; 

This model is integrable and solvable by the same techniques as the hermitian 
one. The conserved invariant quantities are 

I n = trL n = tr(-R) n = tiiiU^Uy 1 (117) 

and the solution is 

U = Be iCt (118) 
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with B a unitary and C a hermitian matrix satisfying 



BCB~ 1 -C = J (119) 
For the Sutherland case with J = £{uv) — 1), Ui — 1, B, C become 

ii 

B jk = 5 jk e^ , C jk = 5 jkPj + (1 - Sjk) ei{q ._ qk) _ 1 (120) 

where, clearly, qi and Pi are initial positions and momenta. 

We conclude by mentioning that, upon scaling x — > ax, t — > a 2 t, the Sutherland 
model goes over to the free Calogero model. This is the 'infinite volume' limit. 



5.2 Unitary matrix model: quantization 

The quantization of the unitary matrix model can be performed in a way practically 
identical to the hermitian model. Indeed, it should be obvious that the matrices 
L and R defined in the two models have idential properties and the whole analysis 
can be done in a parallel way. The new element of the unitary model is that the 
hamiltonian is proportional to the quadratic Casimir I2 = trL 2 , rather than the 
linear one, and appropriate bases of states have to be found in the degenerate spaces 
of 1 1 = trL that diagonalize also I<i- 

We shall, however, give an independent treatment of the unitary model. The 
reasons are primarily pedagogical: we shall use an explicit index notation, rather 
than the neat Tyi calculus of the previous sections, just to familiarize ourselves with 
the alternative. Further, we shall make specific use of the connection of the N x N 
unitary model with the group U(N), taking advantage of (and making contact with) 
known representation facts for this group. 

We begin, again, by defining a canonical momentum matrix conjugate to the 
'coordinate' U 

n = = -u^uu- 1 (121) 

dU 

The Poisson brackets are 

{UjfaUlm} = S jm Sik (122) 

n is somewhat unpleasant, being neither unitary nor hermitian. We prefer to work 
in terms of the hermitian matrices L and R defined previously 

L = iJJU- 1 = -iUU , R = -iU^U = iUU (123) 

Using (j!22)l we derive the following Poisson brackets: 

{Ljk, L im } = i(Lj m Sik — Sj m Li k ) (124) 
{L jk ,R lrn } = (125) 

{Rjk, Rim} = i(Rjm^lk — SjmRlk) (126) 
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The above is nothing but two copies of the U(N) algebra in disguise. To see this, 
expand the matrices L and R in the basis of the fundamental generators of SU (N) 
T a plus the unit matrix: 

N 2 -l 

L = L° + 2 La T a ( 127 ) 

a=l 
N 2 -l 

R = R° + 2 Ra T a ( 128 ) 

a=l 

with L°, L a , R°, R a numbers. Then use the SU(N) commutation relations 

[T a ,T b ] =if abc T c (129) 

as well as the normalization 

ti{T a T b ) = i<U (130) 
to show that the expansion coefficients satisfy the Poisson algebra 

{L a ,L b } = f abc L c (131) 
{L a ,R b } = (132) 
{R a ,R b } = f abc R c (133) 

while L°, R° are central. Note that the U(l) charges 

L° = ~R° = ^{-iU-'U) = ^E*i ( 134 ) 

i 

are essentially the total momentum of the system. 

We are now ready to perform quantization. In the [/-representation, where states 
are functions of U, II becomes the matrix derivative Ujk = —iSu, acting as 

8utr{UB) = B , 5 u tr(U- 1 B) = -U^BU' 1 (135) 

where B is a constant matrix, and similarly on expressions containing more C/'s. L 
and R, upon proper ordering, are represented as 

L = -U5u , R = 6 V ■ U (136) 

where in R it is understood that we first act with the derivative and then right- 
multiply the result by U. With this ordering, L and R become the proper U(N) 
operators acting as 

Ltr{UB) = -UB , LtiiU^B) = BU' 1 (137) 
Rtr(UB) = BU , RtiiU^B) = -U~ l B (138) 
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It is also useful to express their action on arbitrary functions of U as 



tv(ieL)f(U) = /((l - ie)U) - f(U) (139) 
U{ieR)f(U) = f(U(l + ie)) - f(U) (140) 

where e is an arbitrary infinitesimal hermitian matrix, emphasizing their role as 
generators of left- and right-multiplication on U. Correspondingly, the operators L a 
and R a satisfy the SU(N) algebra. Their action can be obtained by taking e = eT a 
with e an infinitesimal scalar parameter, that is, 

isL a f{U) = /((l — ieT a )U) — f(U) (141) 
ieR a f(U) = f{U(l + ieT a ))-f(U) (142) 

The hamiltonian, being classically the kinetic term on the manifold of unitary ma- 
trices U(N), quantum mechanically becomes the laplacian operator on the manifold 
[371] . Using (I127I128|) it is expressed as 

H = hrL 2 = Y^{L a f + l -N{L°f = ^(R a ) 2 + 1 -{R°) 2 = ±trR 2 (143) 

It is, therefore, the common quadratic Casimir of the left- and right-5'L r (iV) algebra 
plus the square of the U(l) charge, the two parts identifiable as the relative and 
center-of-mass energy respectively. 

Quantum mechanical states grouping into irreducible representations of the L 
and R SU(N) algebras will, thus, be degenerate multiplets of the hamiltonian. The 
U(l) (center of mass) part trivially separates: we can boost any state by any desired 
total momentum NP by multiplying the wavefunction by (det U) p . We will examine 
only the SU(N) part from now on. 

A natural basis of states for the Hilbert space are the matrix elements of the 
unitary irreducible representations (irreps) of SU(N). Let R denote such an irrep, 
R(U) the matrix that represents U in this irrep and R a /3{U) the a(3 matrix element 
of this matrix. Clearly a and j3 range from 1 to the dimensionality of R, d R . R a p{U) 
are a complete orthonormal basis of wavefunctions for U, that is 

J [dU\R a p{U)R'^{U)* = 8 RRI 5 ai 5 ps (144) 

with [dU] the volume element on the space of SU(N) matrices as implied by the 
metric ds 2 = —tr(U~ 1 dU) 2 , also called the Haar measure. 

We will, now, show that each R a p(U) is an eigenstate of the hamiltonian with 
eigenvalue equal to the quadratic Casimir of R, Cr. Qualitatively, after the discus- 
sion of the last paragraphs, this should be obvious: L and R generate the trans- 
formations U — > V~ l U and U — > UW. R{U) transforms in the conjugate irrep R 
under L and in the irrep R under R. Since H is the common quadratic Casimir of L 
and R we conclude that all d 2 R states R a p(U) are energy eigenstates with eigenvalue 
Cr = C R . 
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(If you are confused about L generating U — > V~ l U rather than U — > VU, think 
of the difference between active and passive transformations, which is relevant when 
shifting from classical to quantum: ip(x — a) shifts the wavefunction by +a. Also, 
although classical transformations on U compose properly, 

Vx(V 2 U) = (V^U (145) 

quantum mechanically the operators V that perform the shift U — > VU on the 
argument of the wavefunction would compose 

VxiVifiU)) = Vf(V 2 U) = f(V 2 VtU) = {V 2 V x )f{U) (146) 

Therefore we need to invert the action of V to get the right composition law. 

Let us prove the fact HR a p(U) = CrRo^U) more analytically Since R(U) is a 
representation, it obeys the group property 

R a p(UV) = RcrfRyp(V) (147) 

7 

From (|141j) we have 

(1 +ie)L a R af3 (U) = R aP ({\ - ieT a )U) = R ai {\ - ieT a )R lP {U) (148) 
= R a p{U)-ieR a ai R lP {U) (149) 

where R a = R(T a ) is the a-th generator of SU(N) in the R representation. So 

L a R aP (U) = -R^RypiU) (150) 

and 

J2(L a ) 2 R aP (U) = R a a ,R a ,sRsp(U) = Y,(R a )lsRs^(U) (151) 

a a a 

The sum J2 a (R a ) 2 appearing above is the quadratic Casimir in the irrep R and is 
proportional to the identity matrix 5 a s- So, finally, 

HR aP {U) = C R R a(S (U) (152) 

Incidentally, the spectrum spanned by Cr for all R is nothing but the spectrum of 
N free fermions on the circle with the ground state energy and the center-of-mass 
energy subtracted, where the lengths Ri of the rows of the Young tableau of R 
correspond to the "bosonized" fermion momenta 

Pi = Ri-i + 1 (153) 

and where the center-of-mass energy has been subtracted. The condition Ri > 
Ri+i for the rows amounts to the fermionic condition > Pi + \. The spectrum of 
the full matrix model, then, is identical to the free fermion one but with different 
degeneracies. 
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We have, therefore, identified all energy eigenstates of the matrix model. It 
remains to implement the quantum analog of the choice of angular momentum J, 
identify the corresponding reduced quantum model, and pick the subspace of states 
of the full model that belongs to the reduced model. 

J obeys itself the SU(N) algebra (it is traceless, no U(l) charge). As in the 
hermitian model, a choice of value for J amounts to a choice of irrep r for this 
algebra. States within the same irrep are related by unitary transformations of U 
and give the same dynamics; they correspond to a unique physical state. Since 
J = L + R, we see that states transforming under (L, R) in the (R, R) irreps will 
transform in the R x R under J. So, only irreps r that are contained in the direct 
product of two mutually conjugate irreps can be obtained for J. This amounts to 
irreps r with a number of boxes in their Young tableau that is an integer multiple 
of N, just as in the hermitian model case. (To get a feeling of this, consider the 
case N = 2. Then J is an orbital-like realization of the angular momentum through 
derivatives of U and clearly cannot admit spinor representations.) 

We must, therefore, project the d R states in R a p{U) to the subspace of states 
transforming as r under L+R. Call G(R, a; R, /3\r, 7) the Clebsch-Gordan coefficient 
that projects these states to the 7 state of r. Then the relevant states for this model 
become 

*r(U) = E R«p(U)G(R, a; R, /3\r, 7) (154) 

The index 7 labeling the states within r, as we argued before, counts the d r gauge 
copies and does not imply a true degeneracy of states. The degeneracy of the states 
produced by each R is, then, given by the number of times that the irrep r is 
contained in the direct product R x R or, equivalently, the number of times that R 
is contained in R x r. Calling this integer D(R,r; R), we obtain for the spectrum 
and degeneracies: 

E R = C R , D R = D(R, r; R) (155) 

In particular, if Dr = the corresponding energy level is absent from the spectrum. 

Concluding, we mention that an approach which also reproduces the spectrum 
and states of the Sutherland model is two-dimensional Yang-Mills theory on the 
circle |31| 13*2] . This approach is essentially equivalent to the matrix model above 
and we will not be concerned with it. 

5.3 Reduction to spin-particle systems 

So we have derived the spectrum, degeneracy and wavefunctions of the matrix model 
restricted to the sector J = r. Classically these restrictions represented free particles 
(J = 0), Sutherland particles (J = £(vv> — 1)) or something more general. What 
are the corresponding quantum systems? 

To find these, let us reproduce here the expression of the reduced hamiltonian 
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in one of these sectors: 



H = J2^ + \j2^fi^-E. (156) 

This expression remains valid quantum mechanically upon a proper definition (or- 
dering) of the operator K. The only residual quantum effect is a constant term 
E Q that comes from the change of measure from the matrix space to the space of 
eigenvalues. 

Let us expand a bit on this without entering too deeply into the calculations. 
(For details se, e.g., [HH].) The Haar measure in terms of the diagonal and angular 
part of U has the form 

[dU] = A 2 [dV] (157) 
where [dV] is the Haar measure of V and A is the Vandermonde determinant 

A = []2sin^^ (158) 

i<j 

To see this, write the 'line element' — tr([/ _1 <i?7) 2 in terms of V and Xi using ()113j) 
and obtain 



ti{u- l duf = J2 dx l - E 4sin2 ^^(y~ ldv )iA v ~ ldv )ii ( 159 ) 

i,3 



This metric is diagonal in dxi and (V~ 1 dV)ij. The square root of the determinant 
of this metric, which gives the measure (volume element) on the space, is clearly A 2 
times the part coming from V which is the standard Haar measure for V. (We get 
two powers of 4 sin 2 x '~ Xj in the determinant, one from the real and one from the 
imaginary part of (V^dV)^, so the square root of the determinant has one power 
of A 2 .) 

To bring the kinetic Xj-part into a 'flat' form (plain second derivatives in Xi) we 
must multiply the wavefunction with the square root of the relevant measure (com- 
pare with the change from cartesian to spherical coordinates in central potential 
problems). The net result is that the wavefunction \l/ in terms of Xj and V is the 
original wavefunction ip(U) of the matrix model times the Vandermonde determi- 
nant. This, however, also produces an additive constant E D which comes from the 
action of the entire Xj-kinetic operator on A. Noticing that A is nothing but the 
ground state wavefunction of iV free fermions on the circle, we see that E is the 
relevant fermionic ground state energy 

N(N 2 - 1) . , 

Eo = 1 24 ^ (160) 

This is the famous 'fermionization' of the eigenvalues produced by the matrix model 
measure. 
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To determine the proper ordering for K we examine its properties as a gener- 
ator of transformations. Since U = VAV' 1 , and J generates U — > V'UV'~ l = 
(V'V)A(V'V)~ 1 , we see that J generates left-multiplications of the angular part V 
of U. K = V~ X JV, on the other hand, generates ng/ii-multiplications of V, as can 
be see from its form or by explicit calculation through its Poisson brackets. As a 
result, it also obeys the SU(N) algebra. Its proper quantum definition, then, is such 
that it satisfies, as an operator, the SU(N) algebra. It clearly commutes with the 
diagonal part Xi and its momentum pi, since it has no action on it. Its dynamics are 
fully determined by the hamiltonian (jl56j) and its SU(N) commutation relations. 

We can, therefore, in the context of the particle model (|15fi|) . forget where K 
came from and consider it as an independent set of dynamical SU(N) operators. K, 
however, obeys some constraints. The first is that, as is obvious from K = V~ l JV, 
K carries the same irrep r as J. The second is subtler: a right-multiplication of 
V with a diagonal matrix will clearly leave U = VAV" 1 invariant. Therefore, this 
change of V has no counterpart on the 'physical' degrees of freedom of the model 
and is a gauge transformation. As a result, we get the 'Gauss' law' that physical 
states should remain invariant under such transformations. Since K generates right- 
multiplications of V, and Ka (no sum) generates the diagonal ones, we finally obtain 

(no sum) Ka = (on physical states) (161) 

(A more pedestrian but less illuminating way to see it is: J = i[Z7 _1 ,Z7], being a 
commutator, vanishes when sandwiched between the same eigenstate of U. Since K 
is essentially J in the basis where U is diagonal, its diagonal elements vanish.) Note 
that the constraint (|161|) is preserved by the hamiltonian ([156)1 . 

The above fully fixes the reduced model Hilbert space as the product of the 
iV-particle Hilbert space times the <i r -dimensional space of K, with the constraint 
(llfilj) also imposed. The further casting of the model into something with a more 
direct physical interpretation relies upon a convenient realization of K. Any such 
realization will do: simply break the representation of SU(N) that it carries into 
irreps r and read off the spectrum for each r from the results of the previous section. 

We shall implement K in a construction a la Jordan- Wigner, as exposed in the 
hermitian matrix model. Let a m i, a mi , m = l,...q, i = 1, ... N be a set of Nq 
independent bosonic oscillators jHZ|: 

[ami^lj] = 5 mn 5ij (162) 

Then 

K a =J2 OmiTtjOmj ( 163 ) 
m=l 

is a realization of the SU(N) algebra. (T^ are the matrix elements of the funda- 
mental generators T a .) The corresponding matrix elements of K are 




(164) 
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Correspondingly, the coefficient of the Sutherland potential in ()156j) is (for % ^ j) 



K-ijKji ^ ] Q'miQ'ni ^nj^rnj ^ ] ^ J rafi"mi (165) 



We already see that the degrees of freedom of K are redistributed into degrees of 
freedom for each particle in the above. Specifically, defining 



^i,mn ^mi^ni ( ^ ] ^si^si ] &mn (166) 

1 \s=l / 

and comparing with ()164|) we see that the Si are N independent sets of operators 
each satisfying the SU(q) algebra. Before expressing K^Kji in terms of the Si let 
us see what the constraint (jl61j) implies: 

Q t If 

Ku = 22 0, m id m i — — 22 a mk a mk — (167) 
m=l m,k 

J2m,k a lnk a mk commutes with all Kij and all S ijrnn . It is, therefore, a Casimir and can 
be chosen as a fixed integer £N equal to the total number operator of the subspace 
of the oscillator Fock space in which the model lives. The above constraint, then, 
implies 

g 

aLi a mi = t (168) 

m=l 

(We see why we had to choose the total number operator to be a multiple of N: 
the operator in ()168j) above is also a number operator and can have only integer 
eigenvalues.) Using this in ()166|) we can express 



and therefore 



Q'mi^'ni Si mn -\- S mn (169) 

Q 



KijKji — ^2 Si jmn Sj jnm + — Si ■ Sj + (170) 



where Si ■ Si = tv(SiSj) is the S'?7(g)-invariant scalar product of the two SU(q) 
'vectors.' We finally obtain the hamiltonian as 



1 1 2S--S- + e{e+q) 

g-E^+aS ; sin uj (mi 

So it is a Sutherland-like model, in which the particles also carry SU(q) internal 
degrees of freedom ('spins') and the potential contains a pairwise antiferromagnetic 
interaction between the spins. 
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It remains to specify the representation in which the SU(q) spins are and find 
the irreps contained in this realization of K, therefore obtaining the spectrum. A 
realization of the form ()164|) for q — 1 in terms of bosonic oscillators contains all 
totally symmetric irreps of S(N) (that is, the ones with a single row in their Young 
tableau). 1)164)1 is essentially the direct product of q such independent realizations, 
so it contains all direct products of q totally symmetric irreps. This produces all 
irreps with up to q rows in their Young tableau, perhaps more than once each. 
The constraint (J168)) . however, implies that the total number of boxes in the Young 
tableau of these irreps is £N. We recover once more the constraint that we derived 
before based on the origin of component of R x R. 

Similarly, the realization (jlfifij) of Si contains all the totally symmetric irreps of 
SU(q). ()168|) implies that the number of boxes of these irreps is equal to £, so the 
spins Si are each in the Mold symmetric irrep of SU(q). Solving this model amounts 
to decomposing the tensor product of these N spins into irreducible components of 
SU(q). Each such component corresponds to a subspace of the Hilbert space with 
a fixed total spin S. This same irrep, interpreted as an irrep r of SU(N), will be 
the corresponding irrep of K, and also of J, and thus will determine the spectrum 
of this sector through (j!55l) . 

Let us elucidate the above by reproducing the two simplest cases: free particles 
and (spinless) Sutherland particles, comparing with the classical treatment. 

a) Free particles correspond to J = K = 0. So there is no spin and no potential 
and we have non-interacting particles. From (j!55j) we see that all Dr are one, and 
thus the spectrum is the free fermion one, as commented before. The matrix model 
naturally quantizes free particles as fermions. 

b) Spinless Sutherland particles correspond, classically, to J = £(vv< — 1). So 
J is rank one (ignoring the trace). Quantum mechanically this corresponds to the 
irrep r of J having only one row and therefore only one independent Casimir. Since 
q in the realization above corresponds to the number of rows, we must have q = 1. 
Spins, therefore, are absent. The strength of the potential becomes £(£ + 1) where 
£N is the number of boxes in the one row of r. By standard Young tableau rules we 
see that the degeneracy Dr is one if the row lengths of R satisfy 



else it is zero. The spectrum of this model is, then, the same as the spectrum of free 
particles but with the selection rule for their momenta 



We recover the 'minimum distance' selection rule of the Calogero model that led 
to the interpretation as particles with generalized statistics! Only, in this case, the 
statistics parameter £+ 1 is a positive integer. The wavefunctions of the Sutherland 
model can be related to characters of the U(N) group (see, e.g., [Ml). 

A Jordan- Wigner realization of K in terms of fermionic oscillators is also use- 
ful and leads to particles with spins interacting via ferromagnetic Sutherland-type 



R i > Ri+i + £ 



(172) 



Pi > Pi+i + £ + l 



(173) 
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potentials. The hamiltonian becomes [32] 



g = E^- 5 E 43in U-,: I™) 

where now the spins are in the £-fold antisymmetric irrep of SU(q). We will not 
elaborate further and leave the details as an exercise to the reader. 

Clearly there are other, more general realizations, involving a mixture of fermionic 
and bosonic oscillators. Such a construction, e.g., involving a single bosonic and a 
single fermionic oscillator produces the so-called supersymmetric Calogero model 
[3*5] . The model of [3E] also falls into this category. 

We close this section by noting that the unitary matrix model can be thought as 
motion on a space of constant positive curvature U(N)/G, where G is the symmetry 
group that leaves the angular momentum J (with which we reduce the system) 
invariant. For the spinless Sutherland model, we have G = U(N — 1) x U(l), where 
U(N — 1) rotates in the N — l directions normal to the vector v and U(l) is the phase 
arbitrariness of v. The spin systems correspond to more general choices of G. A 
similar construction for spaces of constant negative curvature lead to systems with 
a hyperbolic sine-squared interaction [HI]. Such hyperbolic systems can be easily 
obtained by analytically continuing the variables Xi and pi to the imaginary axis, 
which leaves the hamiltonian real. 



5.4 Concluding remarks about the matrix model approach 

In conclusion, the matrix model has provided us with the following: 

1. An augmentation of the permutation group into the SU(N) group and a 
corresponding possibility to define statistics through the irreps of SU(N). 

2. A realization of generalized scalar statistics but with quantized statistics 
parameter I + 1 in terms of the Calogero model. 

3. A realization of generalized 'non-abelian statistics' in terms of particles with 
internal degrees of freedom interacting through a generalized Calogero-type poten- 
tial. 

4. A systematic way of solving the above models. 
What the matrix models has not provided is 

1. A realization of the Calogero model for fractional values of the interaction 
parameter I. 

2. A realization of spin-Calogero systems with the spins in arbitrary (non- 
symmetric) representations. 

3. A control of the coupling strength of the potential for the spin-Calogero 
models. (Note that the coefficient of Si ■ Sj terms is fixed to ±2 and also the 
constant therm i{l + q) is entirely fixed by the spin representation.) 

There exist generalizations of the above models both at the classical [HS] 13*^] 
and the quantum level jlH]. They all share, however, the limitations spelled out 
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above, especially (3). These restrictions are important, in the quest of more general 
statistics but also from the more practical point of view of solving spin-chain models 
with spins not in the fundamental representation, as we will shortly explain. For 
this reason, yet a different approach will be pursued in the next section, namely, the 
operator approach. 

6 Operator approaches 

The matrix model connection provided us with a powerful tool that not only allowed 
us to generalize the notion of identical particle but also led to the full quantum 
solution of a set of spin-generalized Calogero models. 

As noted, however, in the conclusion of the preceding lecture, the matrix model 
fixes the coefficient of the spin-interaction and scalar interaction terms to ±2 and 
±£(£ ± q) respectively. We cannot choose these coefficients at will. 

We would like to have an approach that defeats this restriction and leads to spin 
models with arbitrary coupling strengths. (This is necessary to attack spin-chain 
systems through the infinite-coupling limit trick to be explained later.) Such an 
approach should also be able to bypass the excursion to matrix models and deal 
more directly with these systems in an algebraic way. This will be achieved with 
the exchange operator formalism [11] . 

6.1 Exchange operator formalism 

Consider the operators My that permute the coordinate degrees of freedom of N 
particles in one dimension which could, in principle, also have internal degrees of 
freedom (M for metathesis, to avoid confusion with momenta pi). They satisfy the 
permutation algebra (symmetric group), in particular 



will be called a one-particle operator (even though it may involve the coordinates 
and momenta of many particles). 

We construct the following one-particle operators jllj : 



Mij = Mr/ = M} 3 = Mji 
[Mij, Mm] =0 if k, I distinct 
MijMjk = MikMij if k distinct 



(175) 
(176) 
(177) 



Any operator Aj on the phase space satisfying 



M i:j A k = 



A k Mij if k distinct 



(178) 
(179) 



= Pi + iW ( x i ~ Xj)Mij =p t + Y iWijMi 



(180) 
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We shall view the tT{ as generalized momenta. To ensure their hermiticity the pre- 
potential W(x) should satisfy 

W(-x) = -W(x)* (181) 
We shall construct the corresponding 'free' hamiltonian from 7Tj 

H = Y, \< (182) 
i z 

In terms of the original Pi this hamiltonian will, in general contain linear terms. To 
ensure that such terms are absent we must further impose 

W(-x) = -W(x) = real (183) 

With the above restriction the hamiltonian H and commutation relations of the 7Tj 
become 

K = E W ijk (M ijk - M jik ) (184) 
k 

# = E ^ + E (Wg + K M a) + £ W ijk M ijk (185) 

i i<j i<j<k 

where we defined the three-body potential and cyclic permutation 

W ijk = WijWjk + WjkWki + W^Wij (186) 
M ijk = M i:j M jk (187) 

To obtain an interesting and tractable model, Wij k , which appears in the commuta- 
tor [7Tj, TTj] and also as a three-body potential, should vanish or at most be a constant. 
This leads to a functional equation for W(x): 

W(x)W(y) - W(x + y) [W(x) + W(y)} = const(= W ijk ) (188) 

We present the solutions: 

a) W ijk = -> W(x) = tjx 

b) W ijk = -I 2 < -> W(a:) = £cotx 

c) = +f > -> W(a:) = fcotha; 
Let's examine each case. 

a) In this case the 7ii become 

if 

vr, = Pl + E — M u (189) 

and satisfy 

[^,^■1=0 (190) 
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The 7Tj commute, so we can consider them as independent momenta. (They are 
sometimes referred to as Dunkl operators [41J.) The hamiltonian reads 



H 



E^+E^ 



(191) 



We obtain a Calogero-like model with exchange interactions. Yet it is nothing 
but a free model in the commuting momenta 71$. Integrability is immediate: the 
permutation-invariant quantities 



E 



Its 



(192) 



obviously commute with each other. If we assume that the particles carry no internal 
degrees of freedom and are bosons or fermions then My = ±1 on physical states. The 
model becomes the standard Calogero model and we have proved its integrability 
in one scoop. (You may be left with a question mark: the hamiltonian and the 
other integrals I n become the standard Calogero ones if My = ±1, so these reduced 
integrals will commute on the bosonic or fermionic subspace; but will they also 
commute on the full Hilbert space? Prove for yourself that this is indeed the case.) 

We can also construct harmonic oscillator operators fTJH2]. The commutators 
between Xj and 7r« are 



Defining 



we can show 



[a i} aj] 



1 1 i + eJ2M l j 



-UMi 



(i * h 



1 
1 

V2 



7Ti — lUJXj 



'jTi + lUJXi) 



(i * h 



-ulM k 
[a], a]} = 



(193) 
(194) 

(195) 
(196) 

(197) 

(198) 
(199) 



This is an extended version of the Heisenberg algebra involving the permutation 
operators. The corresponding oscillator hamiltonian reads 



H = E \i. a \ a i + a i a l) = E \p1 + E \^ x i 



+ E 

i<j 



M l3 ) 



(200) 
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and satisfies 

[H, Oj] = ujcii , [H, a\] = ua\ (201) 

This is the harmonic Calogero model with exchange interactions, which becomes 
again the standard model on bosonic or fermionic subspaces for particles without 
internal degrees of freedom. Since 

H = Y, a U + \Nu + l -lu £ M,j (202) 

i i^j 

we see that on bosonic or fermionic spaces the state annihilated by all a« (if it 
exists) will be the ground state. Solving = we obtain for the ground state 
wavefunction 



i<j 

ipF = \{{sgn{x l0 )\x l0 \- l }e-^A (204) 

i<j 

For £ > the bosonic state is acceptable, while for i < the fermionic one is 
acceptable. In the "wrong" combinations of statistics and sign of £ the ground state 
is not annihilated by the aj, but it is still annihilated by all permutation-invariant 
combinations of the c^. 

From (j201jl we see that we can find the spectrum of this model for fermions 
or bosons by acting on the ground state with all possible permutation-symmetric 
homogeneous polynomials in the a\. A basis bor these is, e.g., 

A n = E(4) n (205) 

i 

So the spectrum is identical to non-interacting fermions or bosons, but with a dif- 
ferent ground state energy. For the 'right' combinations of £ and statistics, where 
()204|) are the correct ground state wavefunctions, the ground state energy is 

N N(N -!),„, , , 

E = -u + ^^ J -\£\u (206) 

which is the correct Calogero result. 
Finally, the quantities 

= = (207) 

i i 

can be shown to commute JT] , and therefore this system is also integrable. It is left 
as an exercise to find the commutation relations of the hi and show that [I n , I m ] = 0. 
b) In the case W(x) = £cotx we have 

7Ti = pi + % cot XijMij (208) 
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K = £(A% - M Jik) (209) 
so the momenta are now coupled. The hamiltonian becomes 

H = E y t+z i^)_ e (W f jl + E M J (210) 

■ z i<3 - bill Xjj y Z . <j<fe J 

We obtain the Sutherland model with exchange interactions plus an extra term. 
On bosonic or fermionic states this becomes an overall constant and we recover the 
standard Sutherland model. Again, since H is by construction positive definite, if a 
state satisfying 71^ = exists it will be the ground state. We obtain 

ipB = II I (211) 

i<j 

4>f = Y[sgn(xij)\ sin Xij\ e (212) 

i<j 

which are acceptable for the same combinations of i and statistics as before. For 
both cases = 1 so 

E. = e* N '- 1 '> (213) 
24 v ' 

is the correct Sutherland model ground state energy. The excited states can again 

be obtained in a (rather complicated) algebraic way Finally, the quantities 

2£ 

= + ' E Mii = ft + e- £ — —-Mtj (214) 

can be shown to have the same commutation relations as the hi defined previously 
for the harmonic system. Therefore, the integrals constructed from them 

4 = £tt? (215) 

i 

commute and the model is integrable. 

c) For W(x) = i cothx we have a similar commutation relation and a hamiltonian 

^I^ + E^ + ^^+EM.,) (-) 

i z i<j smn x ij \ z i<j<k J 

This is the inverse-hyperbolic-sine-square model and supports only scattering states. 
Its integrability can be obtained as for the Sutherland model above, or simply as 
an 'analytic continuation' of that model for imaginary period of space. We will not 
examine it any further. 

In conclusion, an exchange-family of models was introduced, solved and related 
to the standard Calogero models in spaces of definite symmetry. It is remarkable 
that all these proofs work directly, and only, at the quantum domain (there is no 
classical analog of My). 
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6.2 Systems with internal degrees of freedom 

We can easily extend the previous results for particles with internal degrees of free- 
dom. For this, assume that the particles are distinguishable or, equivalently, that 
they carry a number q of (discrete) internal degrees of freedom (species) that can 
be used to (partially) distinguish them. Their states are spanned by \x,a >, where 
a — 1, ... q counts internal states. The total permutation operator T^, then is 

Tij = U.jOij (217) 

where crjj is the operator that permutes the internal states of particles i and j. 

Let us, then, simply take states that are bosonic or fermionic under total particle 
exchange: ly = ±1. On such states 

My = ±<Tij (218) 
and the Calogero and Sutherland exchange model hamiltonians become jUj 

H C = Y. \v\ + E \<f*i + E (219) 

B . =E y + E ,J (220) 



We get the Calogero and Sutherland models with spin-exchange interactions. From 
the completeness relation for the fundamental SU(q) generators T a 

i 2 - 1 i i 

E T af3 T *5 = ^agS-tP ~ TT^aP^S (221) 

we deduce the form of the operators cr^ 

(Tij = 2Si • Sj H — (222) 

where S 1 " acts as T a on the internal states of particle i. So the spin-dependent 
interaction coefficient of the potential in the hamiltonian takes the form EE EJ 

Hj 

=F £ (iSi -SjT^+^j (223) 

We have recovered the ferromagnetic and antiferromagnetic spin model of the pre- 
vious section but with arbitrary coefficient! On the other hand, the spins are nec- 
essarily in the fundamental of SU(q). So we have obtained a generalization of the 
coupling constant with respect to the matrix model but a restriction of the allowed 
spins. 



42 



Note that £ here is an arbitrary parameter, while £ in (jl71j) was the size of the 
symmetric representation of SV For £ = 1 and spins in the fundamental, the matrix 
model and exchange-operator model agree. It is interesting to note that we can go 
from ferromagnetic to antiferromagnetic interactions either by changing the sign of 
£ or by changing the statistics of the particles. 

The solution of the above models can be obtained algebraically. For the spin- 
Sutherland model this is rather complicated and is related to the so-called Yangian 
symmetry jlHJ HH1 EH1 EH- For the spin-Calogero model it is easier jS2J. Let us 
concentrate on the model with interaction £{— 2S{-S^£— -) and define the operators 

4 = £(4) n , (A a J = J2(*l) n S? (224) 

i i 

and their hermitian conjugates. These form a complete set for all permutation- 
symmetric creation and annihilation operators for all species of particles. Yet the 
commutators among themselves and with H do not involve i. They create, therefore, 
the same spectrum of excitations over the ground state as N noninteracting bosons 
or fermions with q species. For £ > the ground state is the bosonic one: 

iJ B = Tl\x^e-^A Xs ({a l }) (225) 

i<j 

where Xs is a totally symmetric state in the <Tj. The set of all Xs forms the N- 
fold symmetric irrep of the total spin S = J2iSi- Therefore the ground state is 
(N + q — 1)\/N\(q — 1)! times degenerate. For £ < the above is not normalizable 
any more. But we remember that we can obtain the same model by starting from 
fermions and the opposite coupling — £ > 0. The ground state, then, is of a fermionic 
type 

= E(-i) p fnO n M-^r j >" (226) 

P \ i } i<j 

where P are total particle permutations and a« are a set of fixed values for the 
indices o\ that determine the state. Clearly the ground state will be obtained for 
the minimal total power of appearing above, and that will happen for a maximally 
different set of values CKj. These states form the n-fold antisymmetric irrep of the 
total spin S, where n = N{modq). The ground state is, thus, q\/n\(q — n)\ times 
degenerate. The above spectra will come handy later. 

6.3 Asymptotic Bethe Ansatz approach 

We already mentioned that there are elaborate algebraic approaches to derive the 
spectrum of the spin-Sutherland model, based on the Yangian symmetry. We will, 
instead, take a lower-key approach which reproduces the same spectra and is phys- 
ically more lucid, although not as rigorous. We will take the ABA route. 

Consider distinguishable particles of the exchange-Calogero type without external 
potential, coming in with asymptotic momenta k{ and scattering off each other. 
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Before scattering, their positions are in some definite ordering determined by the 
ordering of their momenta (it is the inverse of that ordering). 

The key observation is that, after scattering, the particles have simply 'gone 
through' each other with no backscattering The impenetrable 1/x 2 potential 
has become completely penetrable in the presence of the exchange term! You can 
prove this fact by examining the asymptotic properties of a simultaneous eigenstate 
of 7Ti, . . . tcn which is obviously an energy eigenstate: at X{ —>■ ±00 the prepotential 
terms are vanishing and we simply have eigenstates of the individual Pi. Since 
there are no pieces with the values of Pi permuted (coming from backscattering) 
we have complete transmission. (To explicitly see how it works, it is instructive 
to consider the two-body problem, decompose it into symmetric and antisymmetric 
parts, scatter and recombine the parts after scattering. A relative phase of n between 
the two parts is what produces the effect.) 

(Puzzle: what happens with the correspondence principle? With ft back in, the 
interaction coefficient is l[l — HMij). How can a term of order % produce such a 
dramatic effect, particles going through each other, in the H — > limit?) 

So the only effect of the scattering is a phase shift of the wavefunction which, as 
we have said, is the sum of two-body phases 



This is true on an infinite space. On a periodic space we can still use the above 
result, together with the requirement for periodicity for the wavefunction, to derive 
the spectrum. This is the ABA method and is expected to reproduce the correct 
results in the thermodynamic limit of many particles at constant density [53J. It 
gives, in fact, the exact answer for the Sutherland model [2. , so we can expect it to 
work also in the present case. For a space of period 27r the result is 



The left hand side counts the total phase picked up by a particle going round the 
space and scattering off to the other particles in the way. are arbitrary integers, 
ensuring periodicity. There are, however, some constraints on the choice of n, that 
are imposed by continuity from the £ = case: 
-If ki < kj then rij < rij 

-If Hi = rij there is a unique solution, that is, ki < kj and ki > kj represent the 
same state. 

These rules are important to avoid overcounting and to discard spurious solutions. 
With these, the spectrum obtained is the same as the one derived with more rigorous 
methods. For the ordering n\ < . . . njy the solution for ki is 




(227) 




(228) 



7T 




(229) 
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and similarly for other orderings. We see that the ABA momenta ki are the same 
as the pseudomomenta that we have previously defined. 

The bottom line is that the spectrum and degeneracies are the same as those 
of distinguishable particles obeying generalized selection rules for their momentum. 
Still, what fixes the degeneracy of states is the different ways that we can distribute 
the particles to the quantum numbers rii, rather than ki (see, especially, the second 
rule above). A state of N particles with the same rij, for instance, is nondegen- 
erate although they have different ki which would seemingly imply a permutation 
degeneracy. 

For particles with spin the construction above, in combination with the trick of 
the previous subsection of starting with fermions or bosons, produces a spectrum 
with degeneracies the same as those of free particles (the are 'free' quantum 
numbers). As argued before, for ferromagnetic interactions we must choose bosons 
and combine their spins accordingly, while for antiferromagnetic interactions we 
must choose fermions. To spell it out, this means the following: 

1. Choose a set of quantum numbers n*. The ordering is immaterial, since we have 
identical particles, so you can choose n\ < . . . un- 

2. Place your particles on these quantum numbers and put their spins in the ap- 
propriate state. For the ferromagnetic case treat them as bosons: the total spin of 
particles with the same rij transforms in the symmetric tensor product of their spins. 
For the antiferromagnetic case treat them as fermions: the total spin of particles 
with the same rij transforms in the antisymmetric tensor product of their spins; 
clearly up to q can have the same n$ in this case. 

3. Calculate the energy of this state in terms of the ABA momenta (|229|) : E = J2i k 2 - 
It should be obvious that similar rules applied to the spin-Calogero system reproduce 
the spectrum derived in the last subsection. This method can be used to calculate 
both the statistical mechanics (large N) of these systems and the few-body spectra. 

6.4 The freezing trick and spin models 

Now that we have a tractable way of solving spin-Calogero systems with arbitrary 
strength of interaction we can introduce the freezing trick jB5 j and deal with spin 
chain models. 

Consider, first, the previous ferromagnetic or antiferromagnetic spin-Sutherland 
model. Take the limit i — > oo. The potential between the particles goes to infinity, 
so for any finite-energy state the particles will be nearly 'frozen' to their classical 
equilibrium positions. In fact, even the excitation energies around that configuration 
will go to infinity: the ground state energy scales like £ 2 , while the excitations scale 
like N£n + n 2 with n some excitation parameter. So, to leading order in £ the 
spectrum becomes linear and of order t. These excitations correspond, essentially, to 
phonon modes of small oscillations around the equilibrium positions of particles. The 
'stiffness' of oscillations is, of course, proportional to the strength of the potential 
I 2 and the spectrum is proportional to the frequency, of order £. 
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The quantum fluctuations of the particle positions in any state will scale like 
the inverse square root of the oscillator frequency, that is, like l/y/i. But, in the 
hamiltonian, the piece coupling the spins to the kinematical degrees of freedom is 
proportional to 1/ sin 2 ^. In the large-£ limit, thus, this term becomes a constant 
equal to its classical equilibrium value; so, in that limit, spin and kinematical degrees 
of freedom decouple. (Note that the spin part is also of order £.) The hamiltonian 
becomes 

H = H s + £H spm (230) 
with Hs the spinless Sutherland hamiltonian and H spin the spin part 

i<j 2 

where the classical equilibrium positions Xj are equidistant points on the circle: 

% = ^ (232) 

The hamiltonian (|231|) above describes a spin chain consisting of a regular pe- 
riodic lattice of spins in the fundamental of SU(q) coupled through mutual ferro- 
or antiferromagnetic interactions of strength inversely proportional to their chord 
distance. It is the well known SU(q) Haldane-Shastry (HS) model [HUH]- According 
to the above, its spectrum can be found by taking the full spectrum of the corre- 
sponding spin-Sutherland model in the large-£ limit, 'modding out' the spectrum of 
the spinless model and rescaling by a factor Each state will inherit the spin 
representation of its 'parent' spin-Sutherland state. So, both the energy and the 
total spin of the states of the HS model can be determined this way. Commuting 
integrals of this model JU] can also be obtained this way [33]. At the level of the 
partition function at some temperature T we have 

From this, the thermodynamics of the spin chain model can be extracted. 

We will not give the details of this construction here. We urge anyone interested 
to solve this way a few-site (two or three) spin chain, see how it works and deduce 
the 'construction rules' for the spectrum of a general spin chain. Let us simply 
state that the many degeneracies of the spectrum of the HS model (larger that the 
total spin SU(q) symmetry would imply), which is algebraically explained by the 
existence of the Yangian symmetry, can, in this approach, be explained in terms of 
the degeneracies of free particles. (The degeneracies are not identical, due to the 
modding procedure, but related.) 

For the spin-Calogero model a similar limit can be taken, scaling also the external 
oscillator frequency as u — > too to keep the system bound. The classical equilibrium 
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positions of this model are at the roots of the iV-th Hermite polynomial. We ob- 
tain, therefore, a non-regular lattice of spins interacting with a strength inversely 
proportional to the square of their distance j^S]- The spectrum of this model can 
be found quite easily with the above method. Again, we refer to the literature for 
details [52EH1EI!. 

In the continuum limit (N — > oo) the antiferromagnetic version of both the above 
models become c = 1 conformal field theories, the HS containing both chiral sectors 
while the inhomogeneous harmonic one containing just one sector. 

Other models exist and can be solved in this spirit: hierarchical (many-coupling) 
models jSH], supersymmetric models [SHI ED], 'twisted' models [HI] etc. All, however, 
work only for the fundamental representation of some internal group. The big, 
important open problem is to crack a particle system with a higher representation 
for the spins and arbitrary coupling strength. If this is done, through the freezing 
trick we will be able to solve a spin chain with spins in a higher representation. This 
is interesting since we could then see if the antiferromagnetic system for integer 
SU(2) spins develops a mass gap, according to the Haldane conjecture 

7 Folding the Calogero model: new systems by 
reduction 

In previous sections we presented mainly two variants of the Calogero model: the 
rational (Calogero proper) and the trigonometric (Sutherland) one. There are other 
possibilities, in which the two-body potential is further generalized to a Weierstrass 
function (and its hyperbolic reduction) as well as systems with 'reflection' symmetry. 
Spin-Calogero models along these generalizations and involving spin 'twists' are also 
possible. 

Rather than giving an independent analysis of these models, we prefer to present 
how they can be obtained as appropriate reductions of standard (spin) Calogero 
models. The advantage is a conceptual unifiation of the various models, and the 
realization that the 'root' model of particles with inverse-square interaction contains 
all the fundamental information about these systems, everything else stemming from 
it as particular sectors and reductions. 

7.1 Reduction of spinless models 

The main idea is to reduce the Calogero model by some of its discrete symmetries, 
akin to the 'orbifold' construction in field and string theory. 

Consider any hamiltonian dynamical sysstem with some discrete symmetries D. 
Its equations of motion remain invariant under the phase space mapping — > D '{(f), 
where <fi are phase space variables. Then the reduction to the invariant subspace 
= D((f>) is kinematically preserved; that is, the equations of motion do not move 
the system out of this subspace. Therefore, reducing the initial value data to this 
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subspace trivially produces a system as solvable as the original one. The motion 
will be generated by the original hamiltonian on the reduced space. 

The starting point will be the plain vanilla scattering Calogero model 

# = E^ + ^E^ (234) 

i=l Z Z i+j X ij 

where = X{ — Xj. The considered symmetries are: 

• Translation invariance T: Xi — > Xi + a, pi — > p« 

• Parity P: x { -> -x,, p, -> -p* 

• Permutation symmetry M: £j — > iM(i), Pi — ► Pm(i) with M any element of the 
permutation group of N particles. 

Other symmetries will not be useful for our purposes. 

A direct reduction of the system by any of the above symmetries does not produce 
anything nontrivial or sensible: <p = T(0) is possible only in the trivial case a = 0, 
while = P(4>) and = M(4>) requires (some) of the particle coordinates to coincide, 
which is excluded by the infinite two-body potential. We get useful systems only 
when reducing through appropriate products of the above symmetries. These are: 

a) D — PM: We reduce by P and a particular permutation: M{i) = N — i + 1 
(or any other in the same conjugacy class). M is uniquely fixed from the fact that 
D 2 (4>) = M 2 (4>) is a pure permutation and for it not to make any two different 
particle coordinates to coincide we must have M 2 = 1. Further, if for any particle 
M(i) = i then D(xi) = —Xi and the corresponding coordinate is set to zero. So 
M(i) = i can happen for at most one % (so that no two or more particle coordinates 
are put to zero). So M must be a collection of rank two (two-body) permutations, 
and possibly a rank one (trivial) permutation if the number of particles is odd, which 
can always be represented in the form above. 

The constraint 

Xi = -x N - i+ i , pi = -pAr_ i+ i (235) 

effectively reduces the original system into two mirror-images. The reduced hamil- 
tonian is 

" k 2 ^ ( Xl - x 3 ) 2 + 2^(x l + x 3 ) 2 + ^- 2 



xt 



N' 



N 
~2 



9=9[\ + 2[^\) (236) 



where an overall factor of 2 has been discarded and [.], {.} denote integer and 
fractional part, respectively. The second term in the potential is the interaction of 
each particle with the mirror image of each other particle; the third part accounts 
for the interaction of each particle with the mirror image of itself, and with a particle 
fixed at the origin by the constraint (for odd N). 
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Parity symmetry persists in the case where an external harmonic oscillator po- 
tential is added to the system, promoting it to the confining, rather than scattering, 
Calogero model: 

N 1 1 a N 1 

^E^ + ^El + E^ 2 (237) 

i=l Z Z i^j X ij i=l Z 

Reduction by D = PM produces an integrable system similar as above with the 
added harmonic oscillator potential. 

b) D — TM: No finite-rank element of the permutation group will do, since 
repeated application of D would eventually lead to Xi = Xi + ma. We overcome this 
by starting with N' particles and taking the limit N' — > oo. We pick the element: 
M{i) = i+N for some finite N which, for infinite N', is infinite-rank. The constraint: 

x i+N = Xi + a , p i+N = pi (238) 

leads to a system consisting of infinitely many copies of a finite system displaced by 
multiples of a. We can parametrize the particle indices by the pair (i,m), i — 1 . . . N 
and m G Z, where the original index is i + mN. The constraint now reads 



x 



i,m 



Xi + ma , p i>m = pi (239) 



The resulting system is infinite copies of an iV-body system. The reduced hamilto- 
nian is 

m=-oo i=i 2 2 mi „ = _ 00 i j (xij + ma — na) 2 

Z m'=-oo [i=l m=-ooi^j l X U + ma J i m^O J 

In the above summation terms with (i — j,m — n) are omitted since they correspond 
to self-interactions of particles that are excluded from the original Calogero model. 
The summation over mf above accounts for the infinite periodically repeating copies 
of the system and can be dropped. The infinite m-summation accounts for the 
interaction of each particle with the multiple images of each other particle and can 
be performed explicitly. We eventually get 

N 1 1 2 

H = T i -p 2 i +-T i - 9 1 x (241) 
2^a2sin 2 vr^ 1 } 

that is, the Sutherland model. In the above we omitted an irrelevant constant term 

2 

equal to gN-^, arising from the last term in the summand corresponding to particle 
terms with i = j, which account for the interaction of each particle with its own 
infinite images. 
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c) We can formally extend the T symmetry to complex parameter a. As long 
as there is a subset of coordinates in the reduced phase space that remains real 
and generates all other coordinates through use of D, we will have a well-defined 
real subsystem. Applying D = TM for infinitely many particles parametrized by a 
double index i,j G Z, for two complex translations a and b, the constraint is 

Xij = x i+N >j + a = Xi ij+ M' + b , pij = Pi+N'j = Pi,j+M> (242) 

and we end up with a finite system with iV = M'N' particles periodically repeating 
on the complex plane. (Clearly the specific separation of N into N' and M' is 
irrelevant; we could take, e.g., N' — N ', M' — 1.) Similarly to the Sutherland case, 
the hamiltonian becomes infinitely many copies of 

The above sum has a logarithmic ambiguity that is easily regulated by subtracting 
the constant gj (ma + nb) 2 from each term. We end up with 

N> x i 

^ = E^ + 2^ 9?( ^ |a ' b) (244) 

i=l i^j 

that is, s model with an elliptic Weierstrass potential. 

d) D\ = PM together with D 2 = TM: This is a combination of (a) and (b) 
above. We work again with an infinite number of particles. We impose two con- 
straints: 

X-j+l-e — i P—j+l—e = ~Pj 

Xj+N , = Xj + a , p j+N , = pj (245) 

where e = 0, 1 (any other choice of e is equivalent to one of these). Parametrizing 
j = i + mN' by the pair i, m, i — 1 . . . N', m e Z, we have 

%N-i+l-e = a — Xi jTn , £>jv-j + i_ e — — Pi 

Xi, m = Xi + ma , p itm = pi (246) 

so we end up with a finite system of iV = [(N' — e)/2] particles. The reduced 
hamiltonian is infinitely many copies of 

H = V- 2 -V 91X2 | \ y ^ 2 

i=i 2 2 i+j a sm 2 2 i+ . a 2 sm 2 

+ W ^ + 1 -T ^ (247) 
+ 2 Va 2 sin 2 vr^ + 2 V a 2 sin 2 tt^ 1 } 
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with 

"1 (N-e)\ „ f (N 



*'=n2 +8 (— ))< j "=H f - 2 i^)j (248) 

The coordinates can all be taken in the interval (0, a/2) and the particles interact 
with their infinite mirror- images with respect to mirrors placed at x = and x = a/2 
and with particles fixed at x — (if e = 1) and at x — a/2 (if 2{(iV — e)/2} = 1). A 
similar construction can be performed with two complex translations, as in (c) plus 
one parity reversal. We obtain a similar model but with elliptic functions appearing 
instead of inverse sine squares. 

7.2 Integrals of motion via reduction 

The above exhausts the possibilities for spinless particles. Before we proceed to the 
more interesting case of particles with spin, it is instructive to demonstrate how 
the above construction reproduces the conserved integrals of motion of the reduced 
system. We will consider case (b), as the most generic, case (a) being rather trivial. 
The Lax matrix of the original scattering Calogero model, as given in an early 



section, is 



a 

Uj = Pi 8 i:i + {1 - 8a)— (249) 

Xj^j 

Traces of powers of L produce the integrals of motion in involution for the model: 

I k = trL k , k = 1, . . . N (250) 

For the system of case (b), we promote the index i into a pair (i, m) and choose 
a?i,m = %i + f«a as in (|2H9jl . The resulting infinite-dimensional matrix Li m j n can be 
thought of as consisting of infinitely many blocks of size N x N, m,n labeling the 
blocks and i,j the elements of each block: 

a 

Lim,jn Pi^mn^ij (1 ^mn^ij) [ 7 \ (251) 

Xij + (m — n)a 

We observe a block 'translational invariance' of the matrix L in the indices m, n, 
which reflects the invariance of the model under a translation by a. Due to this, we 
can trade the pair m, n for a single index m — n 

Lim,jn = L m — n -ij (252) 

and thus L becomes an infinite collection of N x N matrices L n labeled by n. We 
define the Fourier transform L(a): 

L(a)=Y,e inr7 L n (253) 



in terms of which L„ is 



^ f 1 T L{o)e-™ (254) 

Z7T JO 
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The corresponding integrals of motion Ik are traces of powers of L. Denoting by Tr 
the trace in the infinite-dimensional space labeled by i, m and by tr the trace in the 
N- dimensional space labeled by % alone, we have: 



Ik — TrL — tr(L ni _ n2 • • -L nfe _ ni ) 

ni,...rifc 

(TOj = fii — = tr(L mi • ■ • L mk _ 1 L_ mi — m fc _i) (255) 



Hi mi,...m k _ 1 



The sum over ri\ above produces a trivial infinity. This is due to the summation 
over the infinite copies of the system, just as in the case of the hamiltonian, and 
will be dropped. In terms of the Fourier transformed L(a) the reduced become 
simply 

4 = 7T" r^L(a) k da (256) 

Z7T JO 

It is now a matter of calculating L(a). From the Fourier transform 

y 1 = ™ f or o < a < 2tt (257) 

„ n + x sin na 



we obtain for L(a) 



Pidjj + (1 - ; ~r - £ 6i 



asm7r— 



(258) 



where the last diagonal term linear in 7r — a came from terms with i = j, n ^ in 
L Uj ij. The matrix inside the square bracket, apart from this linear part, is the Lax 
matrix L of the Sutherland model: 

~ t7r£ 

Lij = PiSij + (1 - Sij) — : (259) 

a 

Substituting ()258j) in ()25fij) we note that the exponential factors cancel (due to 
x ili2 + ■ ■ ■ x ih _ lix = 0) and we are left with 

h = — j** iiii-t^^X da = Vtr— -^-—L k - S ^- r{ix-a) s da 
2ti Jo V a J ^ s\(k-s)\ 2vra s Jo V ; 

[*/2] (Tll\ 2n 



where = trL fe are the conserved integrals of the Sutherland model. We obtain 
a linear combination of the intergal Ik and lower integrals of the same parity. The 
appearence of the lower integrals originates from the interaction of each particle with 
its own infinite images. We saw an example of such a term in the constant potential 
term that we omitted from the reduced hamiltonian of case (b). In conclusion, we 
have recovered the integrals of the Sutherland model. 
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7.3 Reduction of spin-Calogero systems 

We extend now these considerations to systems of particles with internal classical 
U (n) degrees of freedom. The corresponding starting spin-Calogero system can be 
obtained by nontrivial reductions of the hermitian matrix model, analogous to the 
ones that led to the spin-Stherland model from the untary matrix model. Alterna- 
tively, we may otain the spin-Calogero model by taking the infinite-volume limit of 
the periodic spin-Sutherland model. We shall also consider the classical version of 
the model, in which the spins become classical phase space variables whose Poisson 
brackets generate the U (n) algebra. 
The hamiltonian of this model reads 

*=i:y+lY, tA ¥ A (26i) 

The Si are a set of independent classical U (n) spins of rank one and length £, that 
is, n x n rank-one hermitian matrices satisfying 

tr(^) 2 = I 2 (262) 

and with Poisson brackets 

{(Si)ab, (Sj)cd} = -iSij [(Si) ad 5 cb - 5 ad (Si) cb ] (263) 

This can be thought of as the classical limit of the spin models obtained in the 
quantum treatment of the matrix model. The spin representation was restricted to 
be totally symmetric, which classically translates to the rank-one condition. 

Just as in the quantum case, such spins can be realized in terms of classical 
oscillators: 

{S t ) ab = A1A\, a,b = l...n (264) 

where (Af , Af) are a set of nN independent classical harmonic oscillator canonical 
pairs with Poisson brackets: 

{AlA b j } = i5 ij 6 ab (265) 

and satisfying the constraint 

^2A$A? = £ for alii (266) 

a 

The above model, in addition to the previous symmetries T, P and M, also 
possesses the symmetry 

• Spin rotations U: Si — > USiU -1 , with U a constant unitary n x n matrix. 
(xi,Pi) remain unchanged. 

Again, reduction by this symmetry alone leads to no interesting system (implying 
either U — 1 or Si — 0). Reduction by PUM or TUM, however, much along the 
lines of the previous PM and TM reductions, produces new and nontrivial results: 



53 



e) D — PUM with P and M as in (a) before, and U a unitary matrix satisfying 
U 2 = 1 (this is necessary since P and M are of rank two). The constraints are 

Xi = -x N - i+1 , pi = -p N - i+ i , Si = USN+i-iU' 1 (267) 

The reduced hamiltonian acquires the form: 



i=i 2 2 (xj — Xj) 2 2 ^ (xj + Xj) 2 



where iV 7 = [iV/2], (7 = 2{A^/2} and S Q = U SqII^ 1 is an extra spin degree of freedom. 

The form of the hamiltonian for the reduced model and its physical interpretation 
simplifies with an appropriate choice of basis for the spins: by using the £7-invariance 
of the full model, we can perform a unitary rotation V to all spins Si — > VSiV~ l . 
This transforms the matrix U appearing in ()269|) into U — > V~ X UV . With an 
appropriate choice of V we can always choose U to be diagonal: U = diag(e l< ^ a ). 
Because of the constraint U 2 = 1, this means that U can take the form U = 
diag(l, . . . 1, — 1, . . . — 1) with n\ (n 2 ) entries equal to 1 (resp. —1). So we see 
that the original U(n) invariance of the model has been broken to U(n\) x U(n 2 ). 
If rii = n 2 there is an additional Z 2 exchange symmetry. 

As in the spinless case (a), we could have started with a spin-Calogero model in 
an external oscillator potential (which shares the same U and P symmetries), and 
obtain a model as above with the extra confining harmonic potential. 

f ) D — TUM with T and M as in (b) before, and U any unitary matrix. The 
constraint on the phase space is 

Xi = x i+N + a , pi = p i+N , S i+N = USiU' 1 (269) 

The system becomes, again, infinite copies of a-translated and [/-rotated systems, 
and the reduced hamiltonian is 

h - L 2* + 2 E L ( + ma)2 ( 27 °) 

1=1 «j m=— 00 V *J 1 / 

In the above we cannot drop terms with i = j any more, since they are now spin- 
dependent rather than constant. Only the term (i = j, m = 0) must be dropped 
from the summation as before. 

Again, the form of the hamiltonian for the reduced model and its physical inter- 
pretation simplifies with an appropriate choice of basis for the spins which makes U 
diagonal: U = diag(e^ a ). The trace in ()270j) then becomes 

n 

ti{S t U m S U- m ) = {S i )ab(S j ) ba e-* m +<* (271) 

a,6=l 
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where <p ab = <fi a — 4>b- The m-summation appearing in (|27()jl gives 

00 ir( QJT m <? JT~ m \ n 00 p -im(j) ab 

E , -^1 j = (S t ) ab (S,) ba E ; * , 2 (272) 



II 



E ^(^) (SiUiSj)^ (273) 



a,6=l 



with the potential V^,(x) being 



m=— 00 



We must distinguish between the cases i 7^ j and i = j. For the case i 7^ j the sum 
can be obtained from the x-derivative of (|257j) : 

V ab (x) = —e~ 1 ^ (-^ - iir<p ab cot — - 7r|0 a6 | ] (275) 
ysm ^ a J 

For the case % = j we must omit the term m = from the summation. We obtain 
an x-independent potential: 

In the above we omitted a constant (a, 6- independent) term equal to ^2 which would 
contribute to the hamiltonian a term proportional to J2i trfS 1 ^) 2 . Due to (|262j) . this is 
an irrelevant constant. With the above, the reduced hamiltonian eventually becomes 

N 1 1 1 
H = E 2^ + 2 S E ^(^') {StU&U + - E E Vfc (277) 

i=l i^j a, 6 i a, 6 

This is a model of particles with U(n) spins interacting through [/ (n) non-invariant 
couplings, due to the presence of the matrix V ah . The original global U (n) invariance 
is, now, broken to the diagonal U(l) n part and only the diagonal components S aa 
of the total spin 

S ab = E(^U (278) 

i 

are conserved. The standard U(n) -invariant spin-Sutherland model is recovered 
upon choosing V ab ~ S a t, in which case the sums over a, b above become a normal 
trace. This is achieved by choosing a =constant, that is, U = e 1 ^. 

The above model is, in fact, the same as the classical model introduced by Blom 
and Langmann |38j . and the author of this review [10], in the particle-spin form in 
which it was recast in [10] : 

H = \ EP? + \ E (E V ab (x^(S t U(S 3 ) ba + ^ 2 n l ) (279) 
+ |EE^(M)tt + ■^Y,Vab(q a qb-S aa S bb ) (280) 

i ab ab 
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To fully see the equivalence, we must observe the following: 

1. In the present construction we expressed the hamiltonian in terms of U(n) 
spins Si. In j^Uj it was, instead, expressed in terms of traceless SU(n) spins Si. By 
(I264j) and (J266)) we have trS^ = £, so the relation between the two is 

Si = Si-- (281) 
n 

2. The expression (|28()j) derived in 0Hj was fully quantum mechanical. It can 
be seen that the term £(£ + n) in (|28()j) classically becomes I 2 (n was a quantum 
correction similar to the shift of the classical angular momentum J 2 to J (J + 1)). 

3. For the rank-one matrices Si we have the relation 

(Si)ab{Si)ba = (Si) aa (Si)bb (282) 

4. In j3U] a set of dynamically conserved charges q a were introduced that can be 
chosen to have any value as long as they sum to zero. 

5. In [U31 the particles were taken to move on the unit circle, that is, a = 2n. 
Doing the above substitutions in (|277jl we see that it becomes practically identical 

to (|280|) . The two expressions differ by constant terms depending on the charges q a 
and the diagonal elements of the total spin S aa . Since both of these quantities are 
constants of the motion, the two models are trivially related. 

g) We can, similarly to (c), extend the above construction to two complex trans- 
lations and corresponding spin rotations. Parametrizing again the infinite number 
of particles with a doublet of indices i,j 6 Z the constraints are 

X i,j = X i+N',j + a > Pi,j = Pi+N',j j Si + N'j = U SijU 

Xi,j+M' + b , pij = Pij+M' , Sij + M> = VSijV 1 (283) 



X 



'■J 



and we end up as before with a finite system with iV = M'N' particles periodically 
repeating on the complex plane. Choosing M' — 1, N' — N, the spin matrices 
Si+ m N,n = Si;m,n are now expressed as 

Si- m ,n = U m V n SiV~ n U- m = V n U m S l U- m V~ n (284) 

The corresponding reduced hamiltonian becomes infinite copies of 

N 1 1 00 i-r(TJ rn V nc t V^ n TJ^ mc i \ 

*-&^ £r£&^ (285) 

Since the two space translations and the corresponding particle permutations 
commute, for consistency the two spin rotations must also commute: 

UVSi^nV^U' 1 = VUSi^U-W- 1 (286) 



56 



which implies 

[U- 1 V- 1 UV,S im , n ]= (287) 
For this to hold for generic Si we must require the 'group commutator' matrix 

u = U^V^UV (288) 

to be proportional to the identity matrix. Clearly u satisfies det(u;) = u n = 1, so 
we obtain 

UV = ojVU , uj = e i2n » (289) 

with v an integer < v < n. U and V then satisfy Weyl's braiding condition which 
characterizes a noncommutative ('quantum') torus jfiH] . 

For v — (u — 1) the matrices C/ and V commute. For uj ^ 1, however, [/ and 
1/ become 'clock' and 'shift' matrices. We shall deal with the two cases separately. 

7.4 Twisted elliptic spin-Calogero models 

In this case U and V commute and can be simultaneously diagonalized. Just as in 
the case of a single U, we choose a basis for the spins that diagonalizes both U and 
V to U = diag(4>i), V = diag(9i). The m, n-sums that appear in ()285|) become 

n oo i m $ ab -i n Q ab 

£ (SiUOS/)* £ 7-— — (290) 

where the term m = n = is omitted if i = j. We obtain again a potential V a b(xij), 
for i ^ j, given by the sum 

oo -im<j> ab -ind ab 

m,^-oo (x + ma + nb) 2 
and a spin self-coupling V^b for i = j, given by 

e -im4> ab -inB ab 

Kb= £ (ma + nb)2 ( 29 2) 

(m,n)^(0,0) l ma + ^ D J 

Note that now, due to the presence of the phase factors, these sums are convergent 
and have no regularization ambiguity. The only ambiguous terms, defined modulo 
an additive constant, are the ones with a = b. We will comment on the impact of 
such a regularization ambiguity in the sequel. 

The potential V a b(x) is a modular function on the complex torus (a, b) with 
quasiperiodicity 

V ab (x + a) = e^Vabix) 

V ab (x + b) = e ie «»V ab (x) (293) 
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It has a double pole at x = 0, with principal part 



V ab (x) = \ + 0(x°) (294) 



x 2 



and no other poles in each cell. These properties uniquely define V ab and allow for 
an expression in terms of theta-functions. We put 



V ab (x) = Ae^^ qV \Y ^ (295) 



where gi^ are the as yet unknown zeros of V ab (x) and the theta-functions appearing 
above have complex period T = b/a. This has the right quasiperiodicity under 
x — > x + a. In order to also have the right quasiperiodicity under x — > x + b, qip 
must satisfy 

qi + <?2 = — {aSab ~ h<Pab) = Qab (296) 
271 

and to have the right behavior around x = we must further have 

A = f*f* , (297) 



+ -) a ^( = ^ (296 
91 ) 



The equations ()29fi|) and ()298j) above determine q\ and ^2, while f!297j) then deter- 
mines A. It may be possible to express q±, q2 in a more explicit form, or to recast 
(129 5j) in a form more symmetric in a,b, by using theta-function identities. Finally, 
the self-coupling V ab can be extracted from V ab ( 

V ab = lim (v ab (x) - ^) (299) 



V x 1 . 

To sum up, we obtain a 'twisted' U(n) non-invariant spin-generalization of the 
elliptic model given by a hamiltonian of the form (j277J) but with the potentials 
appearing now being given by ()295|299j) . The U(n) invariance of the original model 
is, again, broken down to the diagonal abelian sungroup U(l) n due to the dependence 
of the potential on a, b. The [/(n)-invariant spin-Weierstrass model is regained for 
4>ab = Gab = 0, that is, trivial matrices U and V . 

We point out that for 9 ab = <p ab = 0, that is, Q ab = 0, the equations for g 12 
(J296I298J) are satisfied for any qi = — q 2 leading to an apparent arbitrariness. As can 
be seen, however, by applying the addition formula 

9l{x + q)9i(x - q)9 4 (0) 2 = 9 1 (x) 2 9 4 (q) 2 - 9 A {x) 2 9 1 {q) 2 (300) 
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this simply amounts to an arbitrary additive constant to the expression for V a b(x) = 
V(x). This corresponds to the need for regularization for this expression in the 
absence of phases, as explained before. (In the case of the Weierstrass function this 
is fixed by further requiring that the 0(x°) part of the function at x = vanish, 
which picks q = ttT/2 and makes 9^{q) above vanish.) We also point out that we can 
pick any of these values for q\ = —q 2 at the limit Q = by appropriately choosing 
the ratio 4> a b/9 a b as they both go to zero. 

The ambiguity of the terms with a = b can be fixed in the same way: we can 
choose phases (p aa ^ 0, 9 aa ^ 0, evaluate the expressions, and then let (p aa , 9 aa — > 0. 
This will lead to arbitrary additive constants C a , depending on the ratio (paa/9 aa as 
we take them to zero. The same constants, however, will appear in both V aa (x) and 
V aa . Their net contribution to the hamiltonian will be 

^ H =\HH C a (S l )aa(S J ) aa + \ £ £ C a (S t ) aa (S t )aa = \ £ C a (S aa ) 2 (301) 
ijtj a * i a * a 

Since the diagonal components of the total spin S are still constants of the motion, 
due to the residual U(l) n invariance, this amounts to the addition of an overall 
constant, and thus leads to systems that are trivially related. 

The same discussion applies if the angles <p a and 9 a coincide for two or more 
values of a belonging to a subspace of indices J, in which case 4> a b = 9 ab = for 
a, b G I. This will result to a constant additive matrix C a b in the potential for this 
subspace of indices J, leading to an extra contribution to the hamiltonian 

AH= J2 C ab S ab S ba (302) 

a,b£l 

Since (p ab = 9 ab = for a, b E /, however, the corresponding subgroup of U(n) 
remains unbroken, and thus the corresponding components of the total spin S ab 
appearing above are constants of the motion. Once again, the arbitrary terms are 
constant and we essentially obtain a unique system. 

Overall, this is a generalization of the spin- Weierstrass model to one involving In 
phases that break the U(n) invariance and promote the potential to a modular func- 
tion. The potential lives on a complex torus in the coordinates, where translations 
around each nontrivial cycle are accompanied by spin transformations. The model 
obtained in (f) can be though of as the limit of the present model with b — ► oo, in 
which case 9 a become irrelevant. 

7.5 Noncommutative elliptic spin-Calogero model 

In this case 

UV = tuVU , uj = e i2 "" (303) 

with v a nonvanishing integer < v < n. 

To proceed, we must identify the form of U, V. We need the irreducible represen- 
tations of the relation ()303j) . Call k the greatest common divisor of v and n. Then 
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n = km and v = kfj,, for relatively prime m, fi. The irreducible representations for 
U, V are m- dimensional 'clock' and 'shift' matrices. By a global U(n) spin transfor- 
mation we can diagonalize either of U, V. Choosing U diagonal, the general form of 
U and V will be the direct sum of k of the above irreducible representations: 

U = diag{e i * a ,...e i * h -*}®u, V = diag{e i6 °, . . . e**- 1 } ® v (304) 

where <p q ,6 q are arbitrary phases, determining the Casimirs U m and V m , and u,v 
are the m-dimensional clock and shift matrices 

u a /3 = u a 5 a /3 , v a p = 5 a+ i t p (mod m) , a, f3 — 0, . . . m — 1 (305) 

So the acceptable U and V depend on 2k arbitrary parameters. 

To take advantage of the form ()304|) for U, V we partition Si into k 2 blocks of 
dimension m x m each by using the double index notation 

{Si)ab = (Si) P a q p , a = pm + a , b = qm + p (306) 
The U(n) Poisson brackets in this notation are 

{(Sif a % (S s )%} = [(5i)S Stf 5 rq - 5 aS 5 rq (5i)3] (307) 
The m, n-sums that appear in ()285j) then become 



g—im<t>p q —in0 pq ^m(a-/3) 

E (x i? +ma + nb)2 (3 ° 8) 

where the term m = n = is omitted if i = j. 

The above gives a potential interaction between particles i and j in the form 
of a modular function in which depends on the spin components of particles i 
and j. To make the noncommutative character of the spin interaction manifest, we 
perform a change of basis in the spin states and define 

TO = E^ (CT+t> (^)a + ^ ( 309 ) 

This is essentially a discrete Fourier transform in the sum of the a, (3 indices of S^. 
(Note that, for m odd, S^p is actually antiperiodic in the index a if (3 is odd, and 
vice versa. Although we could have defined a properly periodic matrix, we prefer 
this slight inconvenience in order to make the ensuing formulae more symmetric.) 
In fact, it will be convenient to assemble the double indices (a,/3) and (m,n) into 
vectors a and rh. Similarly, we define c = (a, b) and 4> p = (0 P , 9 P ). 
The Poisson brackets of the Si are found from ()307j) 



3x/3 ,Z, S7 .„ _Sxj3 _ 



{(^)l 9 , (Sj)p = iSq uj^ 5 ps (S^ - u-— 6 rq (Si) a+ f 



(310) 
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This is a structure extending the Moyal (star-commutator) algebra, the exponent of 
D being the cross product of the discrete 'momenta' a and (3. For (rs) = (pq), in 
particular, it becomes the torus Fourier transform of the Moyal bracket 

(SiYfi = iS^J - [sx0\ u (3ii) 

where 

[x] u = — r (312) 

UJ2 — UJ 2 

is the u;-deformation of x. This is the so-called trigonometric algebra with periodic 
discrete indices [TH] . 

Finally, by inverting (j309|) and substituting in (|308|) . the potential energy W in 
terms of the Si acquires the form 

w = E E &)% wfm (sis) 

i,j a;p,q 

The above includes two-body interactions, for i ^ j, as well as spin self-couplings, 
for i = j, arising from the interaction of each particle with its own images in different 
cells. The two-body potential WZ 9 (x) is 

i .oxm p i4> vq -fn 

W p a \x) = - E 7 (314) 

a w m^ (x + c- m) 2 v 7 

while the spin self-coupling WZ q is 

W p J = - E i ^^7^4r7- (315) 
m ^ (c • m) 2 v 7 

If the above potentials were independent of the U(n) indices a and p, q, the sum 
over U(n) indices in the potential energy expression ()313|) would simply be a U(n) 
trace and would give the U(n) -invariant coupling between the spins of particles i 
and j multiplying the standard Weierstrass potential of the elliptic Calogero model. 
In the present case, however, the above potential is spin-dependent and breaks U(n) 
invariance, introducing a star-product twist in the indices a and phase shifts 4> p in 
the indices p, q. Generically, the U(n) invariance of the original model is broken 
down to an abelian U(l) k , amounting to the transformation 

(503 - ^ (Si)Z ^ (316) 

If 4> p are equal for k' values of p, the remaining symmetry U(l) k ' is enhanced to 
U(k'), corresponding to mixing the corresponding p-components. 

The case uj = 1, 4> p — reduces to the standard spin-elliptic Calogero- Moser 
model. The case m = 1 (and thus u — 1) reproduces the [/(n)-noninvariant model 
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of the previous subsection. The general case with u ^ 1 is a new classical integrable 
model of the spin-Calogero type with a spin-dependent potential which is a modular 
function of the two-body distance x^. 

The sums appearing in (j314j) and (|315J) could in principle have ambiguities due 
to the logarithmic divergence of the summation over the radial coordinate on the 
complex plane. This is, indeed, the case for the standard Weierstrass function and 
a specific prescription is needed to regularize it. Different prescriptions lead to 
different additive constants in the result. In our case, however, the presence of the 
extra phases renders the sums convergent and there is no regularization ambiguity. 

The potentials can be expressed in terms of theta-functions. W^ q (x) is a modular 
function on the complex torus (a, b) with quasiperiodicity 

Wl\x + a) = e-^ +2 -^ a2 WZ q (x) 

W p d \x + b) = e -*«-^ ai W p d \x) (317) 
It has a double pole at x — 0, with principal part 

WZ q (x) = -^ + 0(x°) (318) 
mx z 

and no other poles in each cell. These properties uniquely define W$ q (x) and allow 
for an expression in terms of theta-functions. We put 

WZ q (x) = A e -*S*" ' } Y — (319) 

°i O) 

where Q12 are the as yet unknown zeros of W$ q (x) and the theta-functions appearing 
above have complex period T = b/a. This has the right quasiperiodicity under 
x — > x + a. In order to also have the right quasiperiodicity under x — > x + b, Qi j2 
must satisfy 

Ql + Q 2 = ^^l + Ha-c (320) 
Z7T m 

and to have the right behavior around x = we must further have 

+ fiM = _ 2t t ai (321) 

A = « ,322) 

ma^ (121)0,(121) 



The equations ()320)1 and (|321j) above determine Qi and Q2, while (132 2j) in turn 
determines A. The self-coupling iy? 9 can then be extracted from W~ q (x) as 



m q = ^AW^)-— 2 ) (323) 
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The sums appearing in (J314)) and ()315j) are in general convergent, due to the 
presence of the phases. For uo = 1, however, the phases are absent and terms with 
p = q have an additive ambiguity due to the need for regularization for the expression 
(I314J) . In the theta- function expression this manifests in the fact that the equations 
for Qi t 2 (I32()|321j) are satisfied for any Q% = —Q2- By applying the addition formula 

9i(x + Q)9 1 (x - Q)6 4 (0) 2 = 6 1 (x) 2 6 4 (Q) 2 - B A {xfBx{Qf (324) 

this is seen indeed to amount to an arbitrary additive constant to the expression 
for W pp (x). The same holds for terms p,q for which <p pq = 0. Such arbitrariness, 
however, corresponds to trivial redefinitions of the model by addition of constants 
of motion, as explained in the previous subsection. 

In conclusion, we identified an integrable generalization of the elliptic spin model 
which breaks the spin U(n) invariance and promotes the potential to a modular 
function introducing noncommutative spin twists. 

8 Epilogue 

This concludes our promenade in Calogero land. We have seen and touched many 
aspects of these models, but have by no means exhausted them. There are various 
other issues that have not been visited. A few of these are listed below. 

There are systems with nearest-neighbor interactions that can also be obtained 
and solved in the exchange operator formalism [S3] . Such systems were not analyzed 
here. We also mention the alternative operator approaches to analyze such models 
used in 

Spin models that can be solved with the techniques outlined in this review include 
many examples not treated here. Apart from the supersymmetric models mentioned 
earlier, we also list the models treated in [HE]. 

The continuous limit of a system of Calogero particles can be described in the 
collective field theory formulation [HZ]. The system exhibits interesting soliton and 
wave solutions |SE1 IMj, whose chiral nature and connection to Benjamin-Davis- 
Acrivos-Ono solitons is an intresting issue [70] . 

Finally, there are certainly lots more topics related to the Calogero model, and 
related references, to which this report has not made justice. The most egregious 
omissions will hopefully be rectified in the next revision. 
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